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Abstract: We present a new reflection process X in a smooth, bounded 
domain D that behaves very much like obliquely reflected Brownian motion, 
except that the direction of reflection depends on an external parameter S, 
which we call spin, that is only allowed to change when the main process is 
on the boundary of D. The model is inspired in a spinning Brownian ball 
that bounces of a moving wall. 

We show strong existence and uniqueness of the joint process (X, S), 
and prove that it has a unique stationary distribution by describing an 
alternative construction of the process (X , S) based on excursions of X 
from dD, and an associated exit system. 
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1. Introduction 

Let D C R™ be a bounded C 2 domain, and let B t be an n-dimensional Brownian 
motion. A pair [X t ,St) with values in D x W is called spinning Brownian 
motion (sBm) if it solves the following stochastic differential equation 



where L t is the local time for X tl and n is the interior unit normal to dD. Our 
assumptions on the coefficients are as follows: 

• cr(-) is an (n x n)-matrix valued, Lipschitz continuous function, and is 
uniformly elliptic, that is, there is a constant c\ > such that £ T cr(:r)£ ^ 
ci |£| 2 for all £ € R", and all x E~D. 

• 7(x, s) = n(x) + t{x, s) is defined for x € dD and s £ W, where ft is the 
interior normal to dD, and f is a Lipschitz vector field on dD x MP such 
that n(x) ■ f(x, s) = for all x € dD and s GM. P , 

• g(-) is a Lipschitz vector field on dD with values in MP. 

• Research conducted at the University of Washington as part of the author's dissertation, 
and partially funded by the NSF grant number DMS 090-6743. 
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dX t 
dS t 



a(X t )dB t +i{X u S t )dL, 
[g(X t )-S t ]dL u 



(1.1) 
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The process X t behaves just like a Brownian diffusion inside D, and is re- 
flected instantaneously in the direction 7 = n + f once it hits the boundary. The 
challenge is that the direction of reflection depends on an external parameter 
St that we call spin, which is updated every time the main process X t hits the 
boundary of D. 

This type of process arises naturally from a physical model that might be 
useful for applications: consider a small ball that spins and moves around a 
planar box following a Brownian path. On the boundary of the box, we put tiny 
wheels which rotate at different speeds, modifying the spin of the ball as well as 
pushing it in a certain (non-normal) direction. In this context, it is natural to 
think of the boundary wheels as an external forcing system that is not affected 
by the hit of the ball: every wheel on the boundary rotates at a speed only 
dependent on its position. 

Even though our inspiration for the model comes from the spinning ball 
bouncing off of a moving boundary, from the mathematical point of view, it is 
natural to regard the process (X, S) as a multidimensional reflected difussion in 
D x R p with a non-elliptic infinitesimal generator. Indeed, setting Z = (X, S) 
we can write (1.1) as 

dZ t = a (Z t )dW t + K(Z t )dL t , (1.2) 

where Wt is an (n + p)-dimensional Brownian motion, <7q(x,s) is the square 
matrix of size n + p obtained from a(x) by augmenting it with zeroes, and 
k(x, s) — (j(x, s),g(x) — s). Since we have d(D x W) = dD x W, the local times 
of (1.1) and (1.2) are the same because Z is in the boundary of its domain if 
and only if X £ dD, and the interior normal to D x W is just (n, P ) where P 
is the zero vector in M. p . Equation (1.2) does not fall within the domain of the 
submartingale problem of Stroock and Varadhan [19] since the diffusion matrix 
a e is not elliptic. It is clear that the trajectories of a given solution to (1.2) 
could correspond to different driving Brownian motions as long as the first n 
components of those Brownian motions coincide, since the trajectories of the 
process Z are not affected at all by those components. In particular, it would 
be impossible to have strong solutions with respect to the filtration generated 
by W in (1.2), and so equation (1.1) better represents our intuition. 

An alternative to the classical submartingale approach was introduced by 
Lions and Sznitman in Theorem 4.1 in [16], where existence of reflected diffusions 
driven by a general semimartingale was shown, but that result only holds for 
smooth, bounded domains. Their approach is based on an analytical solution 
to the deterministic Skorokhod problem (see also [6] and [7, 8] for some non- 
smooth cases), but it does not yield many probabilistic results, such as the Feller 
property, that we need to study spinning Brownian motion and its stationary 
distribution. 

In [2] , a reflected process with inert drift is studied, and existence is obtained 
by constructing a reflected Brownian motion and then adding a drift through a 
Girsanov transformation. One could regard such process as a reflected diffusion 
with non-elliptic generator, but their treatment of existence differs considerably 
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from ours, which we proceed to describe: we taking a mixed approach that 
combines elements from the submartingale and Skorokhod problems to obtain 
strong existence and uniqueness of spinning Brownian motion as a limit of the 
following family of obliquely reflected processes with elliptic generators: 



This pair of equations can be seen to fit within the classical framework of the 
submartingale problem and strong existence and uniqueness follows from stan- 
dard techniques. We name this family of reflected diffusions the approximating 
processes. For each e > 0, the diffusion (X £ , S £ ) is supported in the unbounded 
set Dxl f , and the existence of a stationary distribution for (X e ,S € ) is not 
guaranteed. Nonetheless, we are able to prove that a stationary distribution 
for the approximating process exists through estimation of the moments of S s . 
Uniqueness in this case follows from standard arguments that use the ellipticity 
of the diffusion matrix. 

To show tightness of the family of approximating processes as e — ¥ 0, we 
show how to control the modulus of continuity of the local times associated to 
the family (X £ , S £ ) by using repeatedly Ito's formula for some well chosen test 
functions. Some ideas in this proof have been taken from [19], and have been 
adapted to work with the supremum of increments up to time t. We then are 
able to prove that the limit of the family is a weak solution to (1.1). By using 
some methods that we learned from [16], we prove that all solutions to equation 
(1.1) satisfy pathwise uniqueness, which combined with weak existence yield 
strong existence and uniqueness by well known arguments. 

One of the main objectives of our research is to understand the structure of 
the stationary distribution of spinning Brownian motion. To obtain existence, we 
show that the spin St eventually hits and stays within a certain compact, convex 
set, which is independent of the starting position of the process. A classical 
result for Feller processes then yields the existence of a stationary distribution. 
The most challenging part is to prove that spinning Brownian motion admits 
a unique stationary distribution under the following crucial assumption on the 
vector field g: 

Al There are p + 1 points xi, . . . ,x p +\ on the boundary of D such that for 
every y G W, there exist non negative coefficients Xj such that y — 



We start with an intermediate result that apparently has little to do with the 
stationary distribution, and it is interesting on its own. We identify the compo- 
nents of an exit system (see Section 3.4, or [17] for a definition) for excursions 
away from the boundary, in terms of the local time Lt of the process, and an 
excursion measure that has been constructed in the build up for Theorem 
7.2 in [4] . It has been pointed out to us that it is possible to use the exit system 
(Lt, H x ) to construct a stationary distribution for the process (X, 5 1 ), in a sim- 
ilar manner that is done in Theorem 8.1 in [11]. We do not need to use Doob's 



( 



dXf = <r(X £ )dB t + 7(X|, SfidL.. 
dS e t = edW t + [g(Xf) - Sf] dL t . 



(3.3) 
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machinery to obtain a stationary distribution as spinning Brownian motion hap- 
pens to be a Feller process that eventually stays within a fixed compact set, and 
thus existence of a stationary distribution follows from standard results. 

Our proof of uniqueness of the stationary distribution is an adaptation, and 
somehow a generalization, of an analogous result for Brownian motion with inert 
drift, recently proved by Bass, Burdzy, Chen and Hairer in [2]. Although the 
literature in stationary distributions of elliptic reflecting diffusions is vast, we 
have not found other models or results for process that are similar to spinning 
Brownian motion. Nonetheless, some of the results available for elliptic reflected 
diffusions have been inspiring to understand the challenges of our research. The 
reader can consult the articles [12, 13, 14] to obtain an idea of the treatment of 
the problem in the case of an elliptic generator, and to see hence the differences 
with our approach to the problem. 

One result that helps characterize the stationary distribution was developed 
by Weiss in his unpublished thesis [20] , by a test that involves only the candidate 
to stationary measure, the infinitesimal generator of the diffusion, and the vector 
field defining the directions of reflection. A recent extension of this result to some 
non-smooth domains was carried out by Kang and Ramanan in [15]. Both results 
ask for the submartingale problem associated to the diffusion to be well-posed, 
but they do not request ellipticity of the infinitesimal generator, and thus apply 
to our setting. We make use of this characterization of the stationary distribution 
to produce an explicit example of stationary distribution is a specific case of 
(1.1), and we also suggest simulations to obtain numerical approximations of 
the stationary distribution. 

1.1. Outline 

The paper is organizes as follows: Section 2 contains an explicit example of a 
sBm and its stationary distribution. We also provide some simulation-generated 
graphs of the spin marginal of the stationary distribution for two different 2- 
dimensional sBm. In Section 3 we present the core results, only providing proofs 
that are short and easy to follow. We leave out the technical lemmas and harder 
proofs for Sections 4.1 to 4.3. 

2. Examples 

In his unpublished thesis [20] , Weiss presents a test to characterize any invariant 
measure of the solution to a well-posed submartingale problem. His test only 
works for a smooth setting, and has been recently extended for a large class 
of non-smooth domains by Kang and Ramanan in [15]. Even though we just 
need Weiss' result given our assumptions on the domain D, the fact the result is 
available for more general ones opens a research line that we had not considered 
before. We lay out these results next. 

Set k(x,s) = (j(x,s),g — s). Since St is bounded in the stationary regime 
we can regard the vector k as bounded. It follows that the following theorem 
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from the unpublished dissertation of Weiss [20] can be applied to our setting: 

Let £ = | J2i,j=i a i,j( x ) dx%x 3 + Ya=i ^( x )^7 De a second order differential 
operator, where a^j and &i are bounded, Lipschitz functions. Assume that a 
bounded, Lipschitz vector field k is given on the boundary of a G 2 (R d ) domain 
G, such that re • h(x) > (3 > for x <E dG. Let <f> be a G 2 (R d ) function defining 
the boundary of G. 

Theorem 2.1 (from [20]). Le£ G be compact in Mr and bj and re as before, 
suppose (ai_j(x)) is bounded, continuous, and positive semidefinite satisfying 
V (f){x) T a{x)V (f>(x) > for x in a neighborhood of dG (i.e. the diffusion has 
nonzero random component normal to the boundary). Suppose that fi is a prob- 
ability measure on G with n(dG) — and 

[ £f(x)n(dx) < (2.1) 
Jg 

for all f € G 2 (G) with V/ • k{x) > for x € dG. Suppose that the submartingale 
problem for a, b and k is uniquely solvable starting from any x G G. Then fi is 
an invariant measure of the diffusion. 

This theorem has been successfully used by Harrison, Landau and Shepp 
[12] to give an explicit formula for the stationary distribution /i of obliquely 
reflected Brownian motion in planar domains, in two cases: (a) the domain 
is of class G 2 (C) and bounded, and the reflection coefficient re has a global 
extension to a G 2 (R 2 ) vector field; and (b) the domain is a convex polygon, and 
the reflection coefficient is constant in each face. Their technique to obtain an 
explicit representation is to assume that fi(dx) = p(x)dx and integrate (2.1) by 
parts to obtain a PDE with boundary conditions for p, and solve such equation. 
Our approach to obtain the stationary distribution for some specific cases of 
spinning Brownian motion is based on the same idea. 

2.1. Spinning Brownian motion in a wristband 

Consider the spinning Brownian motion in the strip [—1,1] x R with periodic 
boundary conditions of period 2-7r. This turns the strip into a compact domain 
and our construction from the previous section can be used to define SBM, and 
to prove that there exists only one stationary distribution. In this example we 
compute explicitly such stationary distribution. 

Consider the function g(x,y) = alm(y) — /31{_ij(y) for positive constants 
a, f3, and t(x, y,s) — \xliiy{y), and the associated spinning BM solving the 
equation 

dX v t = dB y t +n(X t )dL t , 
dXf =dB x t +T(X t ,St)dL t , 
dSt =[g{X t )-S t ]dL t . 

Note that the normal depends only on the y-coordinate, and so X\ has the 
distribution of reflected Brownian motion in [—1,1]. In particular, L t depends 
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exclusively on B\. Also, if we identify the points x and x + 2n, the domain 
becomes a compact space and the existence of a unique stationary distribution 
follows with minor and obvious modifications from our theorem. 

It is clear from the equations that the law of (X, S) starting from (x, y; s) 
is the same as the law of (x + X°,S), where (X°,S) starts from (0, y;s). A 
standard argument then shows that the stationary distribution is invariant under 
translations in the x-coordinate. Thus, the stationary distribution of (X x , X v , S) 
can be obtained directly from that of (X V ,S)., that is, in terms of computing 
the stationary distribution, we can completely disregard the component X x . 

Proposition 2.2. The stationary distribution for the process (X x ,X y ,S) is 
given by the positive, integrable function p(x,y\s) = a(s)y + b(s), where 

«(-) = '-fo-fl 6W = - 1 (2.2) 

Proof. It is enough to show that p(y, s) — a(s)y + b(s) is stationary for the 
process (X V ,S) which is a diffusion that solves a well posed submartingale 
problem in the domain G = (—1, 1) x (— j3, a). Set up this way, theorem 2.1 does 
not apply as G is not of class C 2 . Nonethless, since reflection only happens at 
the boundary of (—1, 1) x R, we can find a bounded domain G' of class C 2 such 
that 

G = (-1,1) x (-/3, a) CGCG'C (-1,1) x R, 

and apply the theorem there for the density p' = pic- Another option is to use 
the recent version of Theorem 2.1 for non-smooth domains due to Kang and 
Ramanan [15]. 

Set K,(y,s) — (n(y),g(y) — s), where n(-) is the negative of the sign function. 
The process (X y ,S) that uniquely solves the submartingale process associated 
to d ViV with boundary condition V/ • n(y, s) > for y £ dG' and / G C 2 (G ). 
Since p 1 = outside of G, the following computation is straight forward by 
integration by parts: 

/• pa pi 

/ d yy f{y,s)p'(y,s)dyds = / d yy f(y,s)p(y,s)dyds 
JG' J-pJ-i 

* l 
d v f(y, s)p{y, s) - f(y, s)a(s) 

The boundary condition V/ • iZ(y, s) > for y = ±1 translates into 

[g(y) - s] d s f(y, s) > sgn{y)d y f(y, s) y = ±1. 

As p > 0, we see that d y f(y, s)p(y, s)\ 1 _ 1 < sgn(y) [g(y) - s] p(y, s)d s f(y,s)\ 1 _ 1 . 
Also, direct computation shows that [g(y) — s] p(y, s) — 2s ^ y/(a — s)(/3 + s) 
at y = ±1. Notice that this term vanishes both at s — — (3 and s = a, and 
its partial derivative with respect to s equals to — sgn(y)a(s) at y = ±1. Doing 



ds 

l 
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integration by parts in s, and using the facts above we obtain: 

/ d yy f(y, s)p'(y, s)dyds < 0, 

JG' 

as desired. □ 

The following graph shows the density p(y, s) from different perspectives. 
We have set a = f3 = 1 to simplify the plotting, but the general shape of the 
graphs is maintained. The reader should notice that the stationary density goes 
to infinity both at (I, a) and (—1,-/3). Heuristically, when the process (X y ,S) 
is near (1, a), the change in spin is little since g{X t ) — S t ~ g{x, l) — a = 0. Thus, 
the spin stays around a for a "long" period of time, and thus the occupation 
measure has a lot of weight around (1, a). A similar situation occurs at (—1, —(3). 
This observation does not seem to generalize trivially to higher dimensional 
cases, but it inspires the examples we show in the next section. 




Figure 1. Graph of the density p(x,s) in (2.2). 



2.2. Higher dimensional spin 

The previous example involves one dimensional spin. We proceed to explore 
several examples of two dimensional spin, and their marginal stationary distri- 
butions to illustrate the impact of different vector fields g{x). Our setting is the 
following: consider the strip domain j) = Ex[-l,l]. We will identify any point 
(x, y) in this domain with all x-translations by 2w, that is, {x, y) = (x + 27r, y). 
Let D be the domain obtained from the strip D after this identification of points. 

Let j(x,y;s) = n(y) + t(s)x, where x is the unit vector in the direction x, 
and n(y) is the normal at y = ±1. Despite this is not strictly speaking a domain 
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in R n+P , our proof of existence and uniqueness of the stationary distribution 
apply with very little modifications. Thus, 

dXf = dBf +T(S t )dL t , 

dXf = dB\ + n{Xf)dL t , (2.3) 
dS t = [i {1} {X y t )g{Xf) + l { _ 1} (Xf)g(Xn - S t ] dL t , 

where Bf is a one e dimensional Brownian motion modulo 27r. 

Since X t is on the boundary of D exactly when the component Xf takes the 
values 1 or —1, we have that X\ is a one dimensional reflected Brownian motion 
in the interval [—1, 1]. This process can be obtained independently of St, and 
thus, the local time L t can be constructed independently of the spin St'. 

dX v t = dBf + dL- 1 - dh\, (2.4) 
where L\ is the local time of the process at y = 1, and L^ 1 is the local time at 

y = -i- 

We do not have explicit examples of the stationary distribution /i in this case. 
Instead, we have simulated spinning Brownian motion for different functions g, 
and computed the average occupation time for the spin component St in order 
to estimate its marginal distribution. Precisely, Corollary (3.25) and Theorem 
IV. 9. 3 in [10] imply that the marginal /is of the stationary distribution of (X, S) 
is t 

H S (U)= lim ~ / P(S u EU)du. 

Even though we have not proved an invariance principle for spinning Brownian 
motion, our approach to estimate /is is to discretize time to sample (Bf,Bf) 
at times discrete times to = 0, t& = kS, where S > is fixed and small, and 
determine the increment B tk+1 — B tk as a two dimensional Normal distributed 
random variable with mean zero and variance Sl2- From equation (2.4) we can 
obtain the local times at —1 and at 1, and proceed to compute Xf and St. 

Measure concentrated near a point. A key aspect of the behavior of the 
stationary distribution can be deduce from the differential equation dSt — 
\g(Xt) — St] dL t , just as in the one dimensional case. Say that g(x) — go is an ex- 
tremal point of and that the surface measure of the set Aq = {x : g(x) = go} 
is positive. Intuitively, this implies that X t spends a lot of local time on Ao 
and so St is frequently pushes towards the value go- Further, if St takes a value 
close to go, then the change dSt is very small, since then g(X t ) — St is small. 
The process then is more likely to stay close to such points than to drift away, 
and it is natural to assume that neighborhoods of such points will have a large 
occupation measure. 

Based on this heuristic, wc simulated the spinning Brownian motion (2.3) 
for t(s) = 1 — |s| 2 , and vector g given by g(x, 1) = ^(1,0) and g(x, — 1) = 
i (cos x, sin x) and obtained the following graphs for the marginal stationary 
distribution of the spin S, where it can be seen that the point (1/2, 0) gets most 
of the weight of the measure. 
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Figure 2. Density concentrated near the point (1/2,0). 



Measure concentrated near the axes. The next example represents a situ- 
ation in which most of the density accumulates near the x and y axes. To force 
the spin to spend most of the time near the axes we will define the function g to 
induce only a vertical change at the top of the wristband, and only a horizontal 
change at the bottom of the wristband: 

g(x, 1) = (0, sinir) and g(x, — 1) = (cosx, 0). 

Since most of the excursions of X t from the boundary of the wristband both 
start and end on the same end of the wristband, the spin is rapidly pushed 
towards the axes. The expected graph of the spin marginal of the occupation 
measures, should show high concentration of density around the axes. This is 
exactly what was found through simulations. 




Figure 3. Contour graph of density concentrated near the axes. 
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3. Main results 

3.1. Obliquely reflected Brownian motion 

Definition 3.1. Let D C R™ be a C 2 (R") domain with interior normal denoted 
by n(x) at a point x on its boundary. Given a Lipschitz, bounded map f : dD — > 
R n such that f{x) and n(x) are orthogonal; a Lipschitz, bounded, and strongly 
elliptic matrix a : D 1" x R ra , and a n— dimensional Brownian motion B t ; 
an obliquely reflected Brownian motion is a strong Markov process Z t 
(with respect to the right continuous, complete filtration generated by B t ) with 
continuous sample paths, satisfying 



where L t is an increasing, adapted process, that only increases when Z t is on 
the boundary, that is, dL t — toD{Zt)dL t . To simplify notation, we will often 
use the shorthand j(x) = n(x) + f(x). 

The process just defined behaves as Brownian motion inside the domain D 
and reflects instantaneously at the boundary, in the direction given by 7. The 
process L t is known as local time and can be thought of as the minimal amount 
of push needed to keep the process within D. Whenever Z t hits the boundary of 
D, infinitely many tiny excursions of Z away from the boundary occur, and so 
local time is indeed accumulated, even though Z spends zero (Lebesgue) time 
at the boundary. 

One very successful way of constructing diffusion processes with boundary 
conditions was developed by Stroock and Varadhan [19]. Their submartingale 
problem proved to be a successful extension of their ideas developed to treat 
the well-known martingale problem. The following survey on the submartingale 
problem is based on their original presentation. 

Let D a non-empty, open subset of R n , such that: 

(i) there exists <j) E C b 2 (R n ;R) such that D = ^ _1 (0, 00), and dD = ^{{Q}). 

(ii) \\Vcj)(x)\\ > 1 for all x € dD. 

The following functions will also be given: 

(i) a : [0,oo) x D — >• A/+(R) which is bounded and continuous, 

(ii) b : [0,oo) x D — ► 1" which is bounded and continuous, 

(iii) 7 : [0,oo) x dD R™ which is bounded, continuous, and satisfies that 
(j(t, x),V(f>(x)) > /3 > for t > and x G dD. 

(iv) p : [0, 00) x dD — ¥ [0, 00) which is bounded and continuous. 

Define, for u > and x E D 



dZ t = <j{Z t )dB t + [n{Z t ) + f{Z t )] dL 



't> 



(3.1) 




(3.2) 
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and, for u > and x £ 3D 



1 1 



d 

dx l 



i=l 



We say that a probability measure P on (Q, 5F) solves the submartingale problem 



We say the the submartingale problem is well-posed if it has a unique solution. 

Theorem 3.1 (Theorem 3.1 in [19]). If the matrix a(-) is uniformly elliptic, 
and under the conditions stated above on the coefficients a, b, 7, p, we have that 
for every x £ D there is a solution P, starting from x, to the submartingale 
problem. 

Theorem 3.2 (Theorem 5.4 in [19]). Assume all the previous conditions on 
the coefficients a, 6, 7,/? hold, and also that 7 is Lipschitz. If either p = or 
p is also a locally Lipschitz function which is strictly positive at each point 
of [0,oo) x dD, then the solution to the submartingale problem is unique and 
depends continuously on t and x. 

In general, continuous dependence of the solution to the submartingale prob- 
lem on t and x follows from uniqueness, as observed in the remark following 
Theorem 5.8 in [19]: it can be proved by establishing compactness, taking a 
subsequence that converges, and identifying any limiting process as the unique 
solution corresponding to the limiting starting point. This is known as the Feller 
property of the process. 

Equation (3.1) fits very nicely the framework of the submartingale problem. 
In the case of our interest, the coefficients in the generator (3.2) do not depend 
on time, and we assume the drift term is zero for simplicity. We will drop the 
subindex u from L u due to the independence of time. We have the following 
theorem, which is standard in the literature. 

Theorem 3.3. Assume that D is a domain of class C 2 (R"), not necessarily 
bounded. Obliquely reflected Brownian motion in D exists, as the unique strong 
solution to (3.1). The process is also square integrable when the initial distribu- 
tion of Zq is square integrable. 





is a P-submartingale for any / £ ([0, 00) x E n ) satisfying 

df 

p^-+B t f>0 on [0,oo) X8D. 
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3.2. The approximating process 

In this section, we find a solution to equation (1.2), which we call the approx- 
imating process. As we will see, the approximating process is an obliquely re- 
flected Brownian motion for which the reflection coefficient might become un- 
bounded. However, when it becomes large it roughly points towards the origin, 
which makes it impossible for the process to escape to infinity. This argument 
is formalized in the following theorem. 

Theorem 3.4. Let D and f be as in Theorem 3.3, and let g be a bounded, Lip- 
schitz vector field on dD. Consider n andp dimensional, independent Brownian 
motion processes B t and Wt- For each e > 0, the equation 



has a unique strong solution satisfying: 

(i) X t ED for allt> 0, 

(ii) L t is a boundary local time process that only increases when X t € dD, that 
is, dL t = l 9D (X t )dL t . 

(Hi) The process has infinite lifetime. 

Proof. Fix e > and let ip n e C%°(W, W) be a function satisfying ip n (s) = s for 
|s| < n, and consider the obliquely reflected Brownian motion Z n = (X n ,S n ) 
in Dxl p solving 



Define the sequence of stopping times T n — {t > : IS"™! > n). It is clear that 
for m > n we have that Z^ T = Z™ AT and T m > T n a.s. Let £ = supT„. It is 
clear that a process Z can be defined for t < £ such that ^*at„ = Ztt\T > ana - 
so Z = (X, S) solves (3.3) for t < £. 

The corollary follows if £ = oo a.s. Let <fi be a defining function for the set D, 
that is G cf{W l ) with D = (j)- 1 ((0,oo)), dD = (j)' 1 ({0}), and |V^(ar)| > 1 
on dD. We have that V0 • 7 — V</> • n > 1. It follows from Ito's formula that 



Since (f> is bounded, the quantity 4>(X t AT m ) — 4>(Xo) is bounded by 2 H^H^. Since 
At/) is also bounded, we arrive at 



dX t = a(X t )dB t + j(X t , S t )dL 
dSt = edW t + [g(X t )-S t ]dL t , 



(3.3) 



dX t = o-(X t )dB t +7(X u S t )dL u 
dS t = edW t + [g(X t ) - i> n (S t )} dL t . 



(3.4) 




rtATm 

C^l + t)- / V0(X u )dB u > L tATm , 
Jo 
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where is a positive constant depending only on <f>. By using the exponential 



martingale exp 
a > 



L aS7(j>(X u )dB u — ^2 f Q \S?4>{X U )\ du we get that for any 



MTV. 



E [e QitAT '-] < exp [C a(l + t)] Eexp ^-a y V^(X„)dB u J 

< exp [C a(l + t)] Eexp ^ jf |V0(x u )| 2 du^) 

< exp (C a(l +t) + a 2 C^) , 
and therefore we obtain the estimate, 

E [ e Qi * AT ™] < c a e Ca{1+a)t , (3.5) 

where the constant C does not depend on the starting point of the process. 

We next deduce a very useful representation of St- By Ito's formula, we have 
for e > and t < £ 

d [e Lt S t ] = e Lt S t dL t +ee Lt dW t +e Lt [g(X t ) - S t ] dL t = se Lt dW t +e Lt g(X t )dL t . 
We obtain for t < £, 



-it 



[ e L "dW u + e- Lt [ g(X u )e L "dL u . (3.6) 
Jo Jo 



Under P^ the first term on the right hand side of last equation is bounded, 
and we can uniformly bound the third term as follows: 



-it 



g{X u )e L "dL u 



-it 



gix^de 1 



< 



\9\\oo e 



-it 



f deh < \\g\ 
Jo 



From these uniform bounds, and by using Ito's isometry below, we deduce that 
for large m and some constant C > 0, 



El S-, 



t/\T„ 



< C + 2E 



,-L„ 



t/\T„ 



(AT™ 



e L "dJU u 



< C + 2! 



tAT„ 



= C + 2E/ e 2i "du < C + 2E [ie 2Lt ] , 



(3.7) 



where the constant C is independent of both the starting point of (X, S) and 
of t. On {£ < i} we have |5 tA T m | /ooasm grows to oo and t stays fixed. If 
P(£ < t) > for some t > 0, we have by the monotone convergence theorem 
that lirrim-^oo E 1 5t A T m | = oo, which contradicts (3.7), and so P(£ < t) = 0. 
Since {£ < oo} = |J n {£ < n}, we are done. □ 

Corollary 3.5. Let (X s , S £ ) be a solution to (3.3). For any t>0, 



S t = e- Lt S + ee~ Lt 



e Lu dW u +e 



-it 



g(X u )e L -dL u . (3.. 
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Proof. Since the life time £ = oo almost surely, the claim follows from (3.6). □ 
Fix T > 0. We would like to point out is a very raw bound for the supremum 
of \Sf\, for times < t < T. It comes out from the well known Burkholder- 
Davis-Gundy martingale inequality. 







2 




< 3|S | 2 e- 2it +3e~ 2Lt 


f e L "dW u 


+ 3e- 2Lt | 


v./) 




Jo 







rt - 2 



\g{X u )\de l 



<3\S \ 



t 

e L "dW u 



+ 3||^L(l-e- £ «) 



If E|jSo| < oo, by estimate (3.5) and the Burkholder-Davis-Gundy inequality, 
it follows that 



E 



sup | & I 

0<t<T 



< OO. 



(3.9) 



Before moving on, we would also like to point ou that the process (3.4) 
converges weakly to the approximating process (3.3). To prove this, we first 
need the following lemma. 

Lemma 3.6 (Echeverria [9]). Let D ([0,oo,G]) denote the space of right con- 
tinuous paths on G endowed with the Skorohod topology. Let P, P™, P™> fc be 
probability measures on D, and let {r^} be a sequence of stopping times on T> , 
such that, for all t > 0, 

(i) P (rfc is lower semicontinuous) = 1 

(ii) lim P(r fc < t) = 0. 

k— too 

Suppose that P n = P n,k on 5F Tj , , that for each fixed k, {^ n ' k } n>1 forms a weakly 
compact set in D, and denote Q k any limit of P n,fe . Suppose that Q k = P on 
y Tk . Then P" tend weakly to P. 

Corollary 3.7. Let i\) n e C£°(W,W) be a function satisfying ip n (s) = s for 
\s\ < n, and consider the obliquely refelcted Brownian motion Z n = (X n ,S n ) 
in D x W solving 



As n 
(3.3). 



\ dX t = a{X t )dB t + j(X t ,S t )dL t 
\dS t = sdW t + \g(X t ) - iP n (S t )] dLt 

oo the process above converges weakly to the approximating process 



Proof. Let P be the law of the approximating process as constructed in The- 
orem 3.4. Let P™ be the law of the process Z n = (X n ,S n ) given in the state- 
ment of this corollary, and let P n,fe be the law of Z n stopped at time Tfc = 
inf {t > : \S?\ > k}. 

It is clear form the first part of the proof of Theorem 3.4 that for m > n > k 
we have that Z tATk = Z™ Tk = Z? ATh for all t > 0, and so P"^ = P m < fe . It trivially 
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follows that P™' converges as n — > oo, and that its limits coincide with P on 

Since each of the processes Z™ and Z t are continuous, the sets {7% < i} are 
closed for each k and t, which means that the stopping times are lower 
semicontinuous almost surely. 

The last condition in Lemma 3.6 (lim/ c _>. 00 P(rfe < t) = 0) has been shown in 
the proof of Theorem 3.4. Our claim follows from said lemma. □ 



3.2.1. Second moment estimates 



The following lemma provides a finer estimate for the second moment of St, 
and will be used later to show existence of a stationary distribution for the 
approximating process. 

Lemma 3.8. Fix a compact set M C W containing a nieghborhood of zero. 
There is a constant C > 0, independent of t > 0, such that for x € D and 
s € M the estimate \St\ 2 < C, holds. 

Proof. Under P XjS we see that both the first and third terms from the right 
hand side of equation (3.8) are uniformly bounded. It is enough to bound the 
expectation of Yf, where Y t — e~ Lt J Q e Lu dW u comes from the middle term in 
(3.8). To this end, we will discretize the Ito integral involved. Let tk — tk/N, 
for k = 0, . . . , N, then 



Yf = lim 

N-too 



N-l 

J2 e -^- L W (w tk+1 - w t , 

.k=0 
N-l 

ST^ -1(T,+-T, t . 1 f^ r Tlr \2 



= ^ E e-* L *- L *>+J {W tk+1 -W tk ) + 

N-l N-l 

+ 2 E E e-^- L ^-^- L ^ (W tj+1 - W tj ) (W tk+1 - W tk ) . 
k=a j=k+i 

where the limits are taken in L 1 (P) sense. 

Let 3^ = a((X u ,S u ) ■ u G [s f t]) be the right continuous, completed aug- 
mentation of the sigma algebra of events happening between times s and t. 
Since L t — L s = L t - S ° S , where 9 is the standard shift operator, we have that 
both L t — L s and Wt — W s are 5^-measurable. Also, since Brownian motion is 
time reversible and its increments are independent, 

E,, s (Wt k+1 -W tk \3i k+1 ) =0. 

By conditioning each term of the double sum above on £F| fc i , we readily 
see that its expectation equals zero. Therefore, by independence of Brownian 



M. Duarte E. /Spinning Brownian motion 



16 



increments, 

E*,a {Y f 2 ) 



N-l 



lim E x s 

k=0 
N-l 

lim V" E x s 

J ^fYl ■ * 



fc=0 



-2(L t -L t )r 



-2(i t -i t ) 



(^ fc+1 -^ tfc ) 2 |^ +1 



(tk+1 — tk) ■ 



Conditioning this time on the sigma algebra 3 r t k+1 , by the Markov property, 



N-l 



E, 



Y f 2 ) = lim V E x . 



k=0 
N-l 



E 



Xt k+i ' St k+i 



e -2L t ~t k+1 



(t 



k+1 



tk) 



< 



lim sup supE Xi 



-2L t - 



Define ip(u) = sup^ s ^x,s [e" 2L -]. It follows fr om the calculation above and 
a simple change of variable that 



E x . s (if) < / tp(u)dt 
Jo 



Our main claim would follow if the integral in the last line remains bounded as 
t — y oo . Since <p(u) > 0, this fact follows from next lemma, which concludes our 
proof. □ 

Lemma 3.9. Let <p(u) = supE^ (e~ 2iu ) . Then, 

x,s 

(i) For all u > 0, ip(u) < 1. 

p oo 

(ii) The integral I tp(u)du is finite. 



2q 



Proof. Let x € D. By Fubini's theorem, for any q > 0, 

/•OO 

E x , s [e~ 2L "] = 2 / e- 2h P x , s [L u <h]dh< F x , s {L u <q)(l- e" 2? ) + e 
Jo 

= 1 - (1 - e- 2q ) P x , s (L„ > (?) , 

where in the second to last step, we have split the integral into two: before and 
after h = q, used the fact that L u is increasing, and that F X:S (L u < h) < 1. 

A result by Stroock and Varadhan [19], Corollary 2.3, says that almost surely, 
for any u>0we have L u > 0. Then, for each x £ D we can choose q > such 
that the right hand side above is strictly less than one. Claim (i) would follow 
from the continuity in [x, s) of the solution to the submartingale problem (see 
[19], Theorem 5.4) in case our process lived in a compact set. Nonetheless, since 
7 is bounded, we have that 



{L u <q}Q\ sup \B V \ < dianxD + 

^0<v<u 
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The set on the right hand side is independent of the starting point (x,s), and 
has probability less than one. Thus, claim (i) is true. 

To see (ii), first notice that for t, u > the Markov property implies that 
tp{u + 1) < <p{u)(p{t). Indeed, 

E x>s [e- 2Lt+ "] = E X>S E X)S [e- 2L ^- 2U \3 t ] = E x>s [E XuSt [e~ 2L "] e~ 2U \ 
< E x>s [p{u)e- 2Lt ] < <p(u)<p(t). 

It follows that In is subadditive and so lim t _). 00 ipit) 1 / 1 exists and equals 
inf t>0 ifit) 1 ^. By (ii), such limit is less than or equal to tp(l) = 1 — 26 < 1. It 
follows that for some to > 0, and all t > to we have ^(t) 1 /* < 1 — 5. Choose 
a > such that 1 — 6 = e~ a . For t > to we have <p(t) < e~ at . Since ip(u) < 1, 
it follows that 

ip(t) dt <t + / e~ at dt < oo. 

Jto 



□ 



Theorem 3.10. Fix e > 0. There is a unique stationary distribution for the 
approximating process (X 6 ,S e ). 

Proof. To find a stationary distribution, we will consider the average occupation 
time measures: 

H x t - S {A) = \ I P*, s {{Xt, SI) eA)du = - \ ( P e (u)l A (x, s)du, (3.10) 
1 Jo 1 Jo 

where A is a Borel set of D x R p , and P £ is the strongly continuous semigroup 
associated to the diffusion (X £ ,S £ ). We assume that the starting point (x,s) 
belongs to a fixed, but arbitrary, compact set M. 

We claim that the family {^' s :t>0,x£D,s£ M} is tight. Indeed, for 
the compact set K = D x {s £ M. p : \s\ < f3} we have by Chebyshev's inequality 



1 r* 

tf> s (K c ) = - / P X , S (\S U \> f3)du< sup ■ 

t JO u£[0,t] 



\Su\ 2 



P 2 

Applying Lemma 3.8, we conclude that given e > 0, there is (3 > such that for 
all (x, s) 6 M, and all t > 1 the bound fif' s {K) > 1 — e holds. The claim now 
follows from Prohorov's theorem. 

Recall that the solution to a submartingale problem is a Feller process (see 
comments right after Theorem 5.8 in [19]). Then, Theorem 4.9.3 in [10] applies 
and it follows that any limit point of {fJ-t ' s } ^ s a stationary distribution. 

Finally, we show that such limit is unique. Let /i be any stationary distribu- 
tion. For 5 > 0, define A{6) = {{x,s) eDxMP: dist(x,9L>) > 6, \\s\\ < 
It is clear that A{6) increases to D x MP as 6 i 0, and so by continuity of the 
measure we have that n(A{6)) f 1. Notice that here we use the fact that 

fi{dD xW)= { F X}S (X t E dD) n(dxds) = 0, 
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since the solution to the submartingale problem spends zero Lebesgue time at 
the boundary. 

It follows that we can choose a 6 > such that fj,(A(S)) > h, which we will 
hx from now on. Choose an open set U whose closure is a compact subset of D 
and an open set V in R p such that U x V C A(6). Let T D be the exit time of 
the process X from D. By stationarity and Fubini's theorem, we have for any 
t > 

n(UxV)> [ P x , s (X t eU,S t eV)fi(dxds) 
Ja(s) 

> [ P xfi (X t eU,s + eW t £V,t< T D ) n(dxds). 

JA(8) 

In the last step, we have used that St = s + sWt for t < T D . Let X D be 
the unique process satisfying dXf — a(X[ > )dB t . It is clear that Xf 3 = X t for 
t < T D , and since Wt and B t are independent, 

H(U x V) > [ Pf {X t D e U) P (s + eW t e V)n(dxds) 

JA(S) 

= / Pf (Xf e U) p(Vst; s, y)dyfi(dxds) 

J A(8) JV 

where ¥ D is the law of X D killed at the boundary of D, and p(t;s,y) is the 
transition density for p-dimensional Brownian motion. For (x, s) € A(6) we have 
< S^ 1 . For y G V, we have that ||y|| < 1, and so we have p(y/et;s,y) > c% 
for some positive constant c± that depends only on e, t, and 5. Let m p be the 
Lebesgue measure in R p , and p D be the density of X D killed upon hitting the 
boundary of D. So far we have proved the following lower bound: 

ft(U x V) > Cl m p {V) I [ p D (t;x 1 z)dz^(dxds). 

It is well know that for every t > 0, the function p D (t; x, z) is strictly positive 
when both x and z are away from the boundary. We have that there are positive 
constants c(t) such that p D (t;x,z) > c(t) for x, z 6 A(S). Setting t = 1 in the 
previous calculations yields that for some constant C3 > 0, which only depends 
on £ and <5, we have //([/ x V) > c^m n+p (U x V) for all open sets U,V with 
J7 x V C A(S). Since the process (X, S) does not spend any Lebesgue time at 
the boundary, it follows that to" +p |£) X rp is absolutely continuous with respect 
to any stationary distribution jjt. 

Let n and v be any two stationary distributions. A well known consequence of 
the ergodic decomposition theorem (see [1] for results in ergodic theory) is that 
any two stationary distributions must be mutually singular, which means that 
there is a set K C W l+P such that n(K) = and v(K c ) = 0. From the previous 
paragraph, it follows that m n +P\ DxRP (K) = and m n+p \ Dx m P (K c ) = 0. This 
contradicts the fact that D is an open set, and so there is at most one stationary 
distribution. □ 
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3.3. Convergence to sBm 

We will next show that as e — > 0, there is a unique limit of the approximating 
processes (3.3) and such limit is a solution to our original equation (1.1). Recall 
the previously introduced notation Z e = (Xf, S%), and notice that to solve (3.3) 
is equivalent to solve dZf = (j e (Z t )dB t + K,(Z £ )dL z t for e > 0, with Zf £ D x R d , 
where 

r ln \ ( 7(X, S) 



a e (x,s) 



el, 



k(x, s) 



This is clear since D x M. d = D x M. d , and the unit normal h(x, s) at (x, s) £ 
d(D x K d ) = dD x M d equals the unit normal n(x) at x £ 3D, for any s £ M p . 
This last identification also allows us to readily see that the local time L t of Xf, 
equals that of Zf . 

In this context, one would be tempted to say that spinning Brownian motion 
(1.1) is the process Z®, but this identification would be very imprecise because 
the driving Brownian motion for sBm in (1.1) is n-dimensional, whereas B t 
above is (n+p)-dimensional. Without this difference in the dimension, we would 
not be able to prove strong uniqueness for spinning Brownian motion. 

Theorem 3.11. Let D C R™ be a C 2 domain. Assume the coefficients cr,T,g 
are of class C 2 . uniformly Lipschitz, and bounded. Let B t be a n-dimensional 
Brownian motion in (17,9^). Then, there exist a (strongly) unique, continuous, 
adapted processes (X,S) satisfying (1.1), with X t £ D. 

We proceed to prove several lemmas that will imply together tightness of the 
approximating processes. The following lemma shows the strong equicontinuity 
of the local times. For proofs of the following lemmas, see Section 4.1. 

Lemma 3.12. There exists a constant C = C(D,a, T,g), independent of e, 
such that for any < 8 < 1 



E sup \L\ -L £ f < CS 1/2 . 
\\t-s\<5 J 



(3.11) 



The estimate found in the previous lemma will be used to prove that the 
family (X s , S e ) is tight. Equicontinuity of X £ follows by a very straightforward 
calculation by using once more the Burkholder-Davis-Gundy inequality: 



Xl\ 2 < 



V D {Xl)dXl 



l D (X £ )a(Z £ u )dB u 



l(Z e u )dLl 



Since both a and 7 are uniformly bounded, by Ito's isometry we deduce that 



E 



sup \X\-X\ 

|t-s|<r5 



< C*E 



sup \B t - B s 

\t-s\<6 



\L e t-L\ 



< C8 1 ' 2 , 

(3.12) 
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for < 5 < 1. 

Unfortunately, we didn't obtain the same type of L 2 estimate for S s . The 
main issue is that we haven't found a successful way to express S e in terms 
of (sub)martingales that would allow us to use the Burkholder-Davis-Gundy 

inequality, in order to estimate the supremum of |S' e |g <t<T in terms of E ^JS^I 2 ^ . 

Instead, we will obtain an estimate in L 1 that is still sufficient to prove tightness 
of the approximating process, as we'll see in the following two lemmas. 

Lemma 3.13. Fix T > 0, and (3 G (0,1/4). There exist constants C — 
C(D,a,T,g) and v > 0, independent of e, such that for any < 6 < 1, 



( 



E, 



sup 

\t-r\<6 
0<t,r<T 



<Or(l + N 2 ) 1/2 <r. 



(3.13) 



Theorem 3.14. 

family 9 M = {P5 
£>([0,oo),R n+ f). 



Let P| s denote the law (X e ,S e ) starting from (x,s). Then, the 
I s : < e < 1, x G D, \s\ < M} is tight in the Skorohod space 



Proof. Since all the measures in 7m are laws of processes starting in a compact 
set, all we need to do is control the modulus of continuity. Notice that for A > 
and 8 G (0, 1) 



sup 

,|t-r|«5 



\Zf 



1, 



Zl\ > A < -E 



sup 

\t-r\<S 



\Xf - XI 



By the remark after Lemma 3.12 and Lemma 3.13, we have that the right 
hand side above is bounded above by Ct,m^~ 1 $ 1 ^ 4 , where Ct,m only depends 
on the coefficients of (3.3), and the constants T and M, of course. Therefore, 

for any \,r) > 0, by choosing 5 < {c^\f\ , qj A 1, we have that 



sup F XiS 

E>0 



/ 

sup 

. \t-r\<6 
\0<r,t<T 



\Zl-Zl\>\ \ < v , 



which is the condition we were looking for. □ 
The next lemma will show that any limit of the tight family (Z e ) is a (weak) 
solution of our original equation (1.1). Also, we show that pathwise uniqueness 
for (1.1) holds, which combined with weak existence gives strong existence and 
uniqueness, proving Theorem 3.11. 

Lemma 3.15. Let (Z™ = Z^") be any convergent subsequence of the family 
(indexed by e) of strong solutions to (3.3), as e n —¥ 0, as in Theorem 3.14- Call 
the limit Z = (X 7 S). Then, there exists a local time process L, such that for 
Z = (X,S), we have that (X,S;L) satisfy (1.1). Also, any solution to (1.1) is 
pathwise unique. 
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3-4- Exit system for excursions away from the boundary 

We introduce the notion of Exit System, first developed by Maisonneuve in [17]. 
Let Z be a standard Markov process taking values in a domain E C R" with 
boundary dE. We attach to E a "cemetery" point A outside of E, and we denote 
by 6 the set of functions / : [0, oo) — > M n U {A} that are continuous in some 
interval [0, C) taking values in R™, and are equal to A in [£, oo). 

An exit system for the process Z from dE is a pair (L£,H X ), where L\ is 
a positive additive functional of Z , and {H x } x£ g E is a family of sigma-fmite 
measures on C such that the canonical process is strong Markov on (to,oo) 
under H^. These measures are called excursion laws. 

Excursions of Z from dE will be denoted e or e s , i.e, if s < u and Z S ,Z U 6 dE, 
and Z t dE for t £ (s,u), then e s = {e s (t) — X t + S , t £ [0,u~ s)} and the 
lifetime of such excursion is given by C( e s) = u — s. By convention, e s (i) = A 
for t > (. 

Let at = inf {s > : L* > t} and let I be the set of left endpoints of all 
connected components of (0, oo) \ {t > : Xt £ dE}. The following result is a 
specialized version of the exit system formula. 

Theorem 3.16 (Theorem 1 in [17]). There exists a positive, continuous addi- 
tive functional L* of {X, S) such that, for every x £ E, any positive, bounded, 
predictable process V , and any universally measurable function f : 6 — > [0, oo) 
that vanishes on excursions et identically equal to A, 



E, 



tei 



(fids 





7 




Jo 



V t H Zt (f)dLl 



(3.14) 



Standard notation is used for H x (/) = j e fdH x . 

The exit system formula provides a technical tool to reconstruct the process 
Z t excursion by excursion. A very nice use of the excursion formula allows us 
to "count" excursion with a given property. For instance, let T be the set of 
excursions from dE starting at a time t £ [a, b) and going through an open set 
U C E. We set V t = 1 and / = lr in the exit formula to obtain 















[ Kz u dL* u 




a<u<b 




J a 



The left hand side is the expectation of the number of excursions in F, which 
can be computed using the exit system (dL*,H x ) according to the right hand 
side. 

For spinning Brownian motion in a domain D, note that excursions of (X, S) 
from dD x M. p correspond to excursions of X from dD as S doesn't change 
within any excursion that goes inside D. In view of (3.14), it is enough to 
consider excursion laws {H x } x€dD for an exit system for X. 
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Theorem 3.17. Let P be the law of Brownian motion killed upon exiting D. 
Define 

H^imipf +AS(x) (3.15) 

and let L t be the local time of (X,S), satisfying equation (1.1). Then H x is 
a sigma-finite measure, strongly Markovian with respect to the filtration of the 
driving Brownian motion B t , and (L t , CiH x ) is an exit system from dD x W 
for the process (X, S), for some constant c\ > 0. 

Notice that the exit system formula does not offer a natural way to normalize 
the measures H x . Moreover, if (A t ,H x ) is an exit system and r)(x) defines a 
positive, measurable function in E, then (77 • At,r](x)~ 1 H x ) also defines an exit 
systems that satisfies (3.14). The excursion measures introduced in the previous 
theorem have been used by Burdzy [4] to establish a canonical choice of an exit 
system for reflected Brownian motion in Lipschitz domains. 

For spinning Brownian motion, excursions from dD start exactly at times 
when the local time increases, and thus it is natural that for some positive func- 
tion -q : 3D -4iwe have that (n ■ dL t ,H x ) is an exit system, because excursions 
of sBm don't look any different from those of reflected Brownian motion. By 
reasoning as above, an exit system for sBm should be (dL t ,r]~ 1 (x)tl x ). The 
theorem then proves that r\ = 1. 



3.5. The stationary distribution 



The main goal of the section is to prove existence and uniqueness of the station- 
ary distribution in a slightly simplified version of the spinning Brownian motion 
(1.1). Consider vector fields f : dD x W -> R n and g : dD -> W that are 
Lipschitz continuous and bounded. Assume the following stochastic differential 
equation has a strong Markov solution 

f dX t = dB t + n{X t )dL t + f(X u S t )dL t , , . 

\ dS t =[g(X t )-S t ]dL t . ^' ibj 

One of the advantages of this equation is that it allows us to solve explicitly 
for St in terms of the process X and its local time L. Indeed, set Y t — exp (L t ), 
and apply Ito's formula to the product Y t S t . Since both St and Y t have bounded 
variation a.s., we get that 

S t = Y t - x So + Yf 1 f Y s g(X s ) dL s . (3.17) 
Jo 

Proposition 3.18. The spin process St in (3.16) is the unique process given by 



S t = e~ Lt S + e~ Lt [ g{X u )e L ^dL. 
Jo 



(3.18) 
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Proof. Let S' t be the process defined by (3.18). We will show that S' t satisfies the 
second equation in (3.16), and conclude by pathwise uniqueness that S' t = St- 

As Y t — e Lt is positive for every t > 0, and almost surely finite, the definition 
of S' t yields d(S' t Y t ) = g(X t )dY t , since both Y t and S' t have paths with finite 
variation. By the product rule, and dYt = YtdLt, we get 



One of the issues with the diffusion (X, S) is the lack of a driving Brownian 
motion for the coordinates related to the spin. At an intuitive level, this means 
that the spin S t could be confined to very small regions of the space, regions 
having Hausdorff dimension less than p, and consequently the support of the 
stationary distribution of the process could be singular with respect to Lebesgue 
measure. To make sure this is not the case, we need to impose some condition 
on the infinitesimal change of St, more precisely, on the function g: 
Al There are p+ 1 points xi, . . . , x p +\ on the boundary of D such that for every 
y G W, there exist non negative coefficients Xj such that y = ^j9( x j)- 

From now on, we assume that Al holds, and we fix the points xi, . . . , x p +± 
that realize it. Notice that if y = d( x j)i then Al implies that — y has an 

expansion with non-negative coefficients, and so we have that for every e > 
there are coefficients rjj > such that = J2^=i Vj9( x j)i an d Y^ P j=\ Vj < £ - 

Lemma 3.19. The set U e = jX^ii Vj9( x j) '■ Vj > 0: Vj < e } * s an °P en 

neighborhood of zero for every e > 0. 

Proof. Since U £ = eU\ it suffices to show that U\ is open. The previous dis- 
cussion shows that € U\. By Al, W — U n <£fiU n , and so one of the sets U n 
contains an open set by Baire's category theorem. As U„ = nUi, we deduce that 
U± contains an open set, which we call V. 

Let w €z V. Since we can write w — ^2 r]jg(xj) with Y] rjj < 1, it follows that 
G Ui — w C U2 = 2U\, and we deduce that U± contains an open set around 
zero, and so does every U e . We need to show that a similar property holds at 
every point z G U\ . Let z = Y] \j(f(xj ) £ Ui with = 1 — ^ ■ Then z + Us is a 
neighborhood of z contained in Ui, which shows that U% is open, as we wanted 
to show. □ 

The convex combinations of g will play a fundamental role in the characteriza- 
tion of the support of the stationary distribution of spinning Brownian motion. 
Such set is commonly referred as the convex hull of {g} = {g{x) : x G dD}: 



The spin process St lives eventually in the closure of this set. To see this at an 
intuitive level, for initial condition Sq we have 



Y t dS' t = g(X t )dY t - S' t dY t = [g(X t ) - S' t ] Y t dL 
so S' t satisfies the second equation in (3.16), as we claimed. 



□ 
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The first term on the right hand side goes to zero and the process eventually 
becomes very close to the one starting at zero, which corresponds to the integral 
term above. The latter clearly belongs to the closure of Hg, as it is a positive 
combination of elements in {g} whose coefficients add up to 1 — Yf . But as 
the discussion following Al shows, the point zero can be written as a positive 
combination with coefficients adding up to F t _1 , and so St eventually belongs 
to the closure of Hg. 

Notice that the set Ui/ 2 from Lemma 3.19 is contained in Hg. The following 
proposition says that the support of the stationary distribution of spinning 
Brownian motion lies within the closure of D x Hg. 

Proposition 3.20. Let St be the spin process in (3.16). Then, for arbitrary 
initial conditions (x,s), there is r < oo such that St G Hg for every t > r 
almost surely. 

Proof. By Corollary 3.23 below the hitting time of U\i2 is finite almost surely 
since Ui/ 2 is a neighborhood of zero. Also, U1/2 Q Hg so we conclude that the 
hitting time r of Hg is finite almost surely. 

For t > t we have that S t = Y t 1 F r S' r + Y t 1 f*g{X u )a!Y u . By continuity, 
S T G Hg, and so for some {xj} and {Aj} with Y^j^j = 1 we have by the 
definition of stochastic integral 

00 N 

St = Yf 1 Y T £ AJSCaJ) + ^ E 9(X t] ) (Y tj+1 - Y tj ) , 
3=1 j=i 

where A = {t\, . . . , ijv+i} is a partition of [r, t]. The sum of the coefficients 
above equals 

Y t - 1 Y T + Y t - 1 (Y t -Y T ) = l, 

and as the stochastic integral with respect to bounded variation processes can 
be computed pointwise, we deduce that St G Hg for all t > r. □ 
We prepare to prove the counterpart of the previous proposition. The proof 
consists of four steps. In the first one (Proposition 3.22), we use a support 
theorem and continuity results for the Skorohod map to show that for any given 
point z G D, T > and e > 0, the probability of (Xt,St) to be in a ball 
of radius e around the final point (z, 0) is positive, no matter what the initial 
position is. In the second step, we use the results of Section 3.4 and excursion 
theory to show how the path of X t can be decomposed into several excursions, 
and how spinning Brownian motion up to the first hitting time of a ball U C D 
can be obtained from sBm conditioned on never hitting U, and adding a suitable 
"last excursion" that hits U. This construction is then used in the third step to 
patch together a spinning Brownian motion from several independent spinning 
Brownian motions V/. In the final step, we show how to condition each of the 
y j 's on hitting the boundary of D only at certain places and deduce from this 
that a component of the spin St has a density with respect to Lebesgue measure. 
This procedure is detailed in Theorem 3.24. 
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Lemma 3.21. Let D, r and g be as above. Let r, T > 0, and z £ D. Assume 
that Al holds. Then, for any (x , s ) £Dx W there is w € C([0,T];R") with 
bounded variation such that there is a unique (x, s) £ C([0,T]; DxW) satisfying 
(a;(0),«(0)) = (xo,flo), (s(T),s(T)) = (ar,0), and for t £ [0,T]. 



a;(t) = a?o + + / 7(x(u), s(u))dl(u) 



o 



s(t) = s + / g{x{u)) - s{u) dl(u), 
Jo 

Here, l(-) is a continuous and increasing function, that only increases when 
x{t) £ 3D, that is l(t) = / Q * l 9D (x(u))dl(u). 

Proof. See section 4.3. 

Proposition 3.22. Let D, f, g; r,T > 0, and z € D as in Lemma 3.21. Then, 
for every {xq,Sq) £ D x W there exists p > such that 

V XOtSo ((Xr, 5 T ) € B{z, r) x 5(0, r)) > q, 

where q depends on Tand r, but is independent of (xg, Sq). 

Proof. Let P be the law of standard Brownian motion in R™. By pathwise 
uniqueness, we know that for a.e. u> € supp (P) there is a unique pair (x, s) € 
C([0,T];D x RP), such that for t € [0,T] 

x{t) = xq + u}(t) + / j(x(u), s(u))dl(u), 



ii 



(t) = s + / g(x(u)) - s(u) dl(u), 







where, /(•) is a continuous and increasing, satisfying i(t) = J Q * la£>(x(w))dZ('u), 
that is, it only increases when £ 3D. It is standard to call this function l{-) 
the local time. 

Let ft be the set of continuous to £ D([0,T];R ra ) such that this uniqueness 
hold. We emphasize that P(f2) = 1 and that the ZU constructed in Lemma 3.21 
belongs to il. Define T in f2 by the assignment uj n- (x, s) as above. We claim 
that r is continuous at uJ, where continuity is taken in the sense of uniform 
convergence in compact sets. Indeed, let ujj £ Q be a sequence converging uni- 
formly in [0, T] to uj. Then, by setting z(t) = (x, s)(t) and ((z) = (j(x, s), g(xj), 
we have that the R n+P — valued functions r\j — {ojj,0) converge uniformly to 
rj = (57, 0). By Theorem 3.1 in [6], we have that the unique solutions (xj, sj) to 
the Skorokhod problem with reflecting vector £ in D x R p , and corresponding 
driving function r/j is relatively compact, and any limit is a solution of the corre- 
sponding problem with driving function UJ. By uniqueness (to £ f2) , we deduce 
that (xj,Sj) — > (x,s) in D([0, T], R n+P ). But as all the involved functions are 
continuous, we actually deduce that the latter convergence is uniform in [0,T]. 
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In particular, there is 8 > such that if u G £l<~)Bc(w, S), then the associated 
solution to the Skorohod problem (x, s) G Bq\ i q T ^{(x,'s), r). Thus we have 

P* ,fco ((^t,S t ) g fl(z,r) x B(y,r)) > F xoM ((x,s) G B c ((x,s),r)) 

> F X0 (ue B c (uJ, 5)), 

which is greater than some positive constant q, independent of Xo, by the support 
theorem of Brownian motion. 

□ 

Corollary 3.23. Let r > and t = inf {t > : S t G 5(0, r)}. Then r is /smie 
almost surely. 

Proof. Let JV[a, b] be the event "S 4 is not in B(0, r) for any t G [a, 6]". Then 
P XiS (r < oo) = 1— lim^oo P XiS (^V[0, &]), where the limit is obviously decreasing. 

By Proposition 3.22 we have that P X:S (N [0, T]) < 1 — q, and a standard 
application of the Markov property shows that s (N[0, nT]) < (1 — g)™, which 
yields P^,;, (r < oo) = 1. □ 

We next proceed to introduce some results about the stationary distribution 
of spinning Brownian motion. Our method is very much an adaptation of the 
proof of Theorem 6.1 in [2]. Such argument involves a decomposition of (the law 
of) Xt in several reflecting processes that are somewhat independent of each 
other. To the reader familiar with excursion theory, "independence" is achieved 
by using suitable exit systems. This decomposition allows us to control both the 
local time and the trajectory of the process before hitting a fixed open set U, 
and deduce that no stationary measure can be null in U . 

Theorem 3.24. Let (X 7 S) be spinning Brownian motion solving (1.1) with 
So = 0. Then, for every T > 0, there is an open set U C H^, and c > such 
that for every open B C U it holds that P IiS (St G B) > c m p (B). 

Proof. See section 4.3. 

Corollary 3.25. Spinning Brownian motion has a unique stationary distribu- 
tion, supported in the closure of D x Hg. 

Proof. Fix T > 0, and U given by the previous theorem. Since S does not 
change when X is inside the domain D and X behaves like Brownian motion 
within excursions, we conclude from the Markov property that (Xt+i, Sr+i) 
has a component with a density with respect to (n + d)-dimensional Lebesgue 
measure on a non-empty open subset of E x U. By Proposition 3.20, we can 
assume that U C Hg. 

We can now combine this with the result of Proposition 3.22 using the Markov 
property to see that for some non-empty set U and any starting point (Xq, So) = 
(zo,ka), the process (X ta ,St ) has a positive density with respect to (n + d)- 
dimensional Lebesgue measure on U under P Zo ,k - This property is generally 
referred to as Harris irreducibility of the process (X, S) . 

From Proposition 3.20, we know that any stationary distribution of (X, S) 
has to be supported in the closure of D x Hg, a bounded set, therefore we deduce 
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that {X, S) has at least one stationary distribution from the standard theory 
of Feller processes (see Theorem IV. 9. 3 in [10]). Let [i be one of them. For the 
open set U in the preceding paragraph 

KU) = [ P*,fc \(X to ,S t0 ) G tl] fi(dxdy) > c m d (U) > 0, 

J E 1 J 

which means that any stationary distribution contains U in its support. This 
contradicts Birkhoff's ergodic theorem in case that more than one stationary 
distribution exist, so there is only one stationary distribution. □ 



4. Proofs 

4- 1. Convergence to sBm 

Lemma 3.12 There exists a constant C — C(D,a,f,g), independent of e, such 
that for any < 5 < 1 



E sup \LI-Ltf) KCS 1 / 2 . 

\\t-s\<6 I 



(4.1) 



Proof. As usual, the letter C will represent a constant that may vary from line 
to line, only depending on D, a, /, g. The letter CV will denote a finite constant 
also dependent on T . 

Let G C£(M. n ) be a defining function for the boundary of D, and define 
$ : M™ x W -> R by $(x,s) = <j>(x). Then $ G C b 2 (R" +p ) and it defines the 
boundary of D x W. By using Ito's formula on both <f> and (j> 2 

<t>(Xl) - 4>{Xt) = J V<j ) T a(Xl)dB u + J L u <f>(X*)du + J V<f> • n(X e u )dL e u , 

(4.2) 

0(X|) 2 - ct>{Xlf = 2 f W T a(X t £ )dS u+ f L u <f(Xl)du. (4.3) 

J s J s 

Therefore, by the Burkholder-Davis-Gundy inequality [18, Chapter 4, thm 4.1], 



E sup |<MA^) 2 - cf>(XI) 2 \ 2 < 4E sup 

\l*-s|<<5 / s|<5 

< CE f sup 

V \t-s\<6 



o-(X^)dB l 



cs 2 



UV 1 '<txr{XZ)\ du +C6 2 , 



which is less that CS, as G C 2 (M. n ) and a is bounded. Recall that <f>(x) > 
in D, and that V<p ■ n(x) > 1 for all x G 9-D. Thus, if t > s, we have from (4.2) 
that 



-2cf>{Xl) (4>(X?) - 4>(Xt )) < C 



J Vct) T a(XI)dB l 



+ C(t-s). 
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Notice that (<p(Xf) -tfXf)) 2 = <f>(Xf) a - f(XI) 2 - 2cj>(XI) (0(Xf) - cj>(Xf)). 
Using this decomposition, the previous estimates, Holder inequality, and the 
Burkholder-Davis-Gundy inequality once more, we arrive at 



E sup \4>{X^) — 4>{X e s )\ < Cb^ + CE sup 

y|i-s|<(5 J \\t-s\<6 

< CSi. 



J V(f) T <j{X £ t )dB v 



2\ 1/2 



This estimate helps to show that the local times L e are (uniformly) equicontin- 
uous. Once more, as V$ • k > 1 on the boundary of the domain, by using (4.2) 
we obtain 



1 9 1 



t 2 

• h{Xl)dLl 



< Mxz) - 0(xi)) 2 + 



V<t> T <j{X^ u )dB u 



+ C(t - s) 2 



Since 1 1 — >• is increasing, we readily get from last estimate and Ito's isometry 
that 



E sup m-Ll\ 2 ) <CS^ 2 . 

\\t-s\<8 J 

Notice that the right hand side does not depend on e, or on the starting point 
of the process. □ 
Lemma 3.13 Fix T > 0, and (3 € (0,1/4). There exist constants v > and 
C = C(D,a,T,g), independent of e, such that for any < S < 1, 



E, 



sup \S\-Sl\ <Ct(1 + M 2 ) 1/2 <T. 



|t-r|<« 
\0<t,r<T 



(4.4) 



Proof. By the product rule we see that d[e Lt Sf\ = ee L tdW t + g(X^)de L t . For 
< r < t < T, we have that 



SI - S s r = ( e -W-i?) - l) S° r + ee- L > J* e L ^dW u + e~ L ' j* 



g{Xl)e L -dL: 



We'll handle the three terms above on the right hand side independently. 

First term: Set E E r t = | (e~ (i ?" i " ) - l) S£|. For x > 0, the function e~ x is 
decreasing, and its derivative is bounded above by 1, thus, for x > y > we 
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have < 1 — e~^ x ~ y ^ < x — y. We estimate the expectation of E$ as follows 

sup E' rt 
\t-r\<5 



< Ex.s 


sup 

\\t-r\<5 


m 




sup \S 

re[0,T] 














1/2 




1/2 


<E X , S 


sup 

\t-r\<6 




-L%\ 2 




sup |^| 2 

re[0,T] 





< cs 



^( N V(l + e 2 )C T ) 1/2 , 



where the bound on the increment of local time comes from Lemma 3.12 and 
the bound on S £ comes from equation (3.9). 

Second term: Let's call this term M£ t , ignoring the e that multiplies it. We 
compute, 



m-K) / e K-K dw 



e L *dW„ o , 



where 9 r is the usual shift operator. By the Markov property, and Doob's mar- 
tingale inequality, 



sup I M e r 1 1 

\t-r\<6 



< E„ 



E Xr ,S r SU P 



e L ^dW u 



< CE Xt 

< CS 1/2 supE^ 



1/2 



,2Lf 



1/2 



The third line follows from Burkholder-Davis inequality, and thus, from (3.5) 
we get the estimate 



E, 



Sup \Mr t \ 
|t-r|<<5 



< C8 1 ' 2 , 5<l. 



Third term: This is the simplest one to handle, by using the fact that the 
local time is increasing and g is bounded 



g{Xl)e L ^dIf u 



<C\L%-L% 



The main claim of this lemma then follows by Lemma 3.12, Holder inequality, 
and the previous estimates for E% t and M^ t . □ 
Lemma 3.15 Let (Z™ = Z^ n ) be any convergent subsequence of the family 
(indexed by e) of strong solutions to (3.3), as e n — > 0, as in Theorem 3.14- Call 
the limit Z = (X,S). Then, there exists a local time process L, such that for 
Z = (X,S), we have that (X,S;L) satisfy (1.1). Also, any solution to (1.1) is 
pathwise unique. 
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Proof. We first prove that any limit satisfies (1.1). By Ito's formula for Zf, we 
have that for any h G C${R n+p ) 



h{Z s t ) - h(Zl) = / l DxKP VM^)a E (^)dB u 



, RP L u h(Zl)du 



Vh ■ it(Z*)dL s u . 



(4.5) 



Our idea is to pass to the limit as e n — > in this equation. First we take care 
of the term involving local time. 

We can actually get stronger convergence of the family of local times Lf 
if we use (4.1) together with Arzela-Ascoli's theorem in the Banach space 
L 2 (C([0, T])). It follows that there is a Cauchy subsequence (£f n ), and L G 
L 2 (C([0,T])), such that 



E sup \L\ n -L s \ z 0. 
\se[o,T] J 

This sequence has a further subsequence that converges pointwise to L. Thus 
the map t L t is increasing and continuous, and also satisfies an estimate 
analogous to (4.1). 

Let / G Cb(R™ +p ) be an arbitrary function. We have 



E 



' l f(Z2)dL2- f f{Z u )dL u 



f(Zl) [dLZ - dL 



< 21 



<2t||/CE sup \Ll-L u \* ) +2E 
\we[o,t] 




[f(zi) - f{z u )] du 



[f(ZD - f(Z u )} dL u 



We know that the first term above goes to zero as n — > oo. The second one 
also does, as f(Z™) — > f(Z u ) pointwise, and L u is continuous with bounded 
variation (see [18, Chapter IV, Thm 2.12]). We use this convergence to prove 
two different things: first, if a sequence of smooth functions tfk 1d ; then, as 
1d(X u ) is locally bounded 







<p k {XZ)dL n u / Vk (X u )dL u , 
o Jo 



so the latter term is a.s. zero. By taking k — > oo, monotone convergence ensures 
that dL t = lg]j(X t )dL t , and so L t is a local time process for X t (and for Z t ). 
Assume that Z% = (x,s) a.s. Then Z = (x,s), and for any h G C$(R n+p ) 
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we have that, almost surely, 

h(Z?) - h{Z$) -> h(Z t ) - h(Z ), 
Vh(ZZ) T a en (Z r :) -> V x h(Z u )a(X u ), 

n-\-d n 

a^Z^didjhiZu) -> ^(^^'(Xujaj^^Zu). 

Also, by setting /(z) = Vh(z) ■ k(z) we have that, almost surely, J V/i • 
K{Zl)dL\ ->■ Jo Vft • n(Z u )dL u . 

It follows from (4.5), and the dominated convergence theorem that the fol- 
lowing equation holds 

h(Z t )-h(Z ) = lim / lz> X RpV/i(3> e „(X£)di? u + 

+ / l_DxRp^°/i(^u)du + / Vh ■ k(Z u )cLL u . 
Jo Jo 

By Ito's isometry, we have 



E 



f* 2 f* 

/ [Vho- en (ZJ)-VxM^»)]^« < / E|V/w 6 „ro- V x ^a(Z„)| 2 du. 
Jo Jo 



Since Wha £n (Z^) — 'Vxha(Z u ) —> almost surely, and ft. € Cg, by the dominated 
convergence theorem in P<g> [0,f] we deduce that the right hand side in the last 
inequality goes to zero. Therefore, 

h{Z t ) - h(Z ) = f l DxMP V x h(Z u )<r(X u )dB u + f l DxRP L a u h(Z u )du + 
Jo Jo 

Vh-K(Z u )dL u , (4.6) 



where L°h — \ ^™,- = i(f T cr)*' J 9i9 7 /i. This proves weak existence for (1.1). 
Proof of pathwise uniqueness. To prove pathwise uniqueness, we will use some 
ideas from [16]. According to Lemma 4.1 in [16], there exists a symmetric-matrix 
valued function A(z), with z £ R n+P , such that A(-) is uniformly elliptic and 
bounded, and A(z) T 7(z) = n(z) for z G dD x MP. Moreover, Ay(-) G C h 2 (K" +p ), 
and for some Co > 

C \z~y\ 2 + K(z) T A{z)(z-y) > z G dD x R p , y G D x K p . (4.7) 

Notice that, for some a to be determined, as k and A are Lipschitz, and bounded, 
we get for z G dD x W, y G D x E p , 

a |z - y| 2 + n{z) T K(y){z ~y)>a\z-y\ 2 - [Cg iA + C ] \z - y\ 2 , 
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thus we get for some C' Q > 0, that C \z — y\ 2 + K(y) T A(y)(z — y) > 0. Adding 
this last inequality with (4.7) together, and putting A(z, y) — A(z) + A(y), we 
finally obtain 

C \z-y\ 2 + K(z)A{z,y){z-y) >0 z e 3D x W ,y eD x K p . (4.8) 

Let 6 be a C 2 function defining dD. Suppose that (X,S;L) and (X',S';L') 
solve (1.1), and set 

r)(z, z') d = exp [\(6(x) + 6{x'))] , 
£(*, z') d =l f n(z, z')(z - z') T A(z, z'){z - z'). 
Let A be a negative number whose value will be specified later. Our aim is to get 



a quadratic bound of the form £(Zt,Z' t ) < M t + C\ Z t 



Z' t \ 



To this end we will 



write the Ito formula for t;(Z t ,Z' t ) and combine terms appropriately to obtain 
terms of order 0(\Z t — Z' t \ ). We start by bounding the Ito terms associated to 
the quadratic variation, for which we introduce the notation z = (x, s), and Ai. 
to represent the i-th row of the matrix A. Similar notation will be used for the 
matrix A. Straight forward computations yields 

d 2 i 



dzJdz 



d 2 tj 
dz'idz 1 



2rjA i3 {z, z) 



= 0(\z-z'\ 
+ 2rj(z,z' 

= 0(\z- z'\')-2 V A lJ (z,z') + 

r fh 

+ 2r)(z,z' 

Interchanging z and z' we obtain aMSm and 



, 86 . . 36 . 



^ 96 A X ^ A 

3z l dz' 1 



3A,. dAi 



dzi dz l 



dA,. dA, 



dz'i 3z % 



(z-z% 



(z-zJ). 



d 2 £ 



The quadratic variation 



3zidz n dz'^dz" 

terms involved in the Ito formula for (,{Z tl Z' t ) are the ones associated to the 
components X t and X' t as the local time is continuous and of bounded variation, 
and they are standard to compute. 
Define, for z, z' E R n+P 

T\ j (z,z') = [a{x)a{x) T ] l3 ^—r- + [a(x)a(x') T ]ij g^jg^ , 



3 2 £ 



d 2 z 



r«(z,z>) = laix'Mx'Y}^^- + [aix'MxYh^-. 

Notice how T2 is obtained from Ti by reversing the roles of z and z' , including 
the variables in the derivatives. Using our previous computations, and using the 
fact that cr(-) is Lipschitz and bounded, we obtain the following equation: 

if (z, z 1 ) =0{\z- z'f) + 2i 1 A lJ (z, z') \a(x) (a(x) - a(x')f 



, , , 86 . dA^ 
2 V (z,z')(X^-A. + — 



(z-z 1 ) <j{x){a{x) +a{x')) J 
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We can obtain an estimate for T2 in a similar way. Adding the two together, 
and using that A{z, z') = A{z' ', z), we arrive at 

if + lf (z,z') = 0(\z-z'f) + 2r 1 A ij (z,z')(a(x) - a(x'))(a(x) - a(x')) T ij + 
+ 2rj(z, z>) (a J^. + ^ (z - z')a{x) (a(x) + a{x')) T l3 + 

- 2 V (z, z>) (x-j^jAi. + (z - z')a(x') (a(x) + a(x')f ij . 

As cr(-) is Lipschitz, the second term above is 0(\z — z'\ 2 ). As both </>(•) and 
A(-) have bounded second derivatives, we have that y X^-(y)A(z, z')i. + 
is Lipschitz, so the second and third line put together are 0(\z — z'\ 2 ) as well. 
Summing up, there is a constant C > 0, independent of t, such that 

i{Z u Z' t )<\ V z Z(Z s ,Z' a )dZ a + f VA{Zs,Z' a )dZ' a + C f \Z S -Z'fds. 
Jo Jo Jo 

Next we take care of the term involving the gradient of £. We have 

H = \^t{z, z') + V (z, z') (2A{z, z% (zj -z' 3 ) + {z- ^') T |^) [z - z'). 

(4.9) 

Since any solution of (1.1) is square integrable, it follows from the boundedness 
of <j>, 77, and A and equation (4.9) that the integral M t — L V Z £,(Z S , Z' s )dB s 
defines a martingale. Using (4.9), the fact that • 7 > 1 on dD, and (4.8) we 
get for A < 0, 

t(Z t , Z' t ) < M t + C [ \Z S - Z'f ds + X f i{Z s ,Z' s )d{L s + L' s ) + 
Jo Jo 

+ 2C [ t](Z„Z' B )\Z s -Z' s \ 2 d{L s + L' s ) 
Jo 

<M t + C I \Z S -Z'fds + 
Jo 

2(vX + C ) / r)(Z s ,Z' s )\Z s -Z'/d(L s + L' s ), 
Jo 

since A is elliptic of constant v > 0, and A < 0. At this point we pick A so that 
vX + Co = 0. Using ellipticity of A once again, we arrive at 

ft 

E\Z t -Z' t \ 2 <C E\Z S - Z'f ds, 
Jo 

and we conclude by Gronwall's inequality that, almost surely, Z t — Z[ for all 
t>0. □ 
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4-2. Exit system for excursions away from the boundary 

Theorem 3.17. Let V D be the law of Brownian motion killed upon exiting D. 
Define 



and let L t be the local time of (X,S), satisfying equation (1.1). Then H x is 
a sigma-finite measure, strongly Markovian with respect to the filtration of the 
driving Brownian motion Bt, and (Lt, ciH x ) is an exit system from 3D x W 
for the process (X, S), for some c\ > independent of x. 

Proof. The fact that H x is sigma-finite for all x £ 3D, and strongly Markovian 
is proved in Theorem 7.2 in [4]. The proof actually applies as these are properties 
of the measures and do not have anything to do with the local time. 

Let (dL't , H x ) be an exit system for (X, S), where L* is the additive functional 
from Theorem 3.16. We will prove that it is possible to replace L\ by the local 
time L t from equation (1.1). Let K C dD be open in the relative topology, and 
set C = H7 — f?|| . Let e > be a very small number. For T > 0, denote by A e 
the set of excursions from 3D that start before time T and reach a level e from 
dD: 

A e d = je t : t < T, sup dist (e t (s), dD) > e j . 

Let x £ 3D be fixed. Pick 5 > so small such that the following is pos- 
sible: Choose a small r and coordinates (x J )™ =1 that differ from the canonical 
coordinates in R n by a translation and a rotation, such that x = in these 
coordinates, and 

• D fl B r (x) — {(y, s) : s > 4>{y)} fl B r (x), that is, D is locally is the graph 
above a function cj>, 

. 3DnB r (x) = {(y,^(y))}nB r (x), 

• V<t>{x) = 0. 

We call the above a normal set of coordinates centered at x. 

Next we decompose K into sets closed set K\, . . . , K ms such that the surface 
measure of the symmetric difference between Kj and Ki is zero (i 7^ j), and 
such that in each Kj we have sets of normal coordinates centered at x\, . . . , x ms , 
and respective boundary defining functions (j) 1 , . . . , 4> ms as above, and 

n(z) ■ n(x 3 ) > I - d 2 y Z £K 3 . (4.10) 

To simplify notation, for z — (y,(j> l ('y)) £ 3D, we will write (^(z) instead of 
(ft (y). This is obviously an abuse of notation that can be justified by identifying 
(fp with (p' ] o i„_i on the boundary of D, where i n -i is the projection into the 
first n — 1 coordinates. By our choice of Vc/^Xj) = 0, we have that nixj) = e n . 
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1 — S < n(z) ■ e„ 



e„ = \/ 1 



1 



We deduce that for small <5, the estimate UV^fz)]! < 28 holds, and by the 
mean value theorem 

\4?(z)\ = \<fP{z) - 4r>(xj)\ < 2S\\xj - z\\ < 2Sr s . 

Let T x . be the tangent plane to dD at Xj . The computation above says that 
the piece of the boundary dD n B rs (xj ) is contained in the cylinder [T x . n 
B rs (xj)] x [—2Srs, 2Srs}. Thus, Kj C Cj D dD. In our proof, it will be useful to 

set A = A(<5) = f 2r$\f8 and define the cylinder Cj — [T Xj n B rs (xj)] x [—A, A]. 
The "top" of the cylinder Cj will be denoted by Cj . 



2 A = ir s VS 




Figure 4. Cylinder Cj, and other related objects. 



By the exit system formula with V u = 1, and the construction above, we have 
that 

E x .s I l K {X u )H Xu (A e )dL* u =E XyS V t K {X u )t A e{e u ) 

Jo u<T 



(4.11) 



u<T 



for any Borel set R C R p , and any T > 0. 



Let cr J ' = inf {t > : X t £ Kj} and T j = inf {t > : X t <£ Cj}. Set a° a = a 3 , 



and for each integer k > 1, let Tj? 



<ri and a{ +1 



<r n o i 



where 6 is the usual shift operator. It is a well known fact that all these objects 
are stopping times. Since Kj is closed, the process lKj{Xu) is predictable as X 
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is continuous. Then, V u — V u lKj {X u ) is non- negative, predictable and bounded, 
and so, by the exit formula for (<£LJ,H X ), and a simple change of variable, 

E x . s V lAj {X u )l A e (e„) = E x<s [ \ Kj {X u )U Xu (A £ )dL* u 

u<T J ° 

GO „ T i 

= ^E IiS l {ffi<T} / + l K] (X u )H Xu (A")dL* u 

oo / T i \ 

= E^V^t} / lK 3 (^)Hx u (A e )dL; o fl(aj) 



OO / „ T 3 



where the last equality holds by the strong Markov property applied at time ai. 

Assume that for all z — (x, s) with x € Cj n d-D, the following equation holds 
for some ci > independent of x 



E 



, / l Kj {X u )U Xu (A s )dL* u = (1 + 0(5))E, / lx 3 .(X u )H Xu (A E )c 1 dL, 1 
Jo Jo 

+ 0(S)E z (L T j), (4.12) 



where O(S) is standard notation for a bounded function that converges to zero 
as S — >• 0. Then, we can trace back all of our computations up to (4.11) to obtain 
the following property: for all measurable sets K C dD and v < t, 

E x . s I l K {X u )H Xu {A e )dL* u = 0(S)E x , s L t + 

+ (1 + 0(5))-E Xi3 f l K (X u ) Cl U Xu (A e )dL u . 

Since K, e > 0, and v < t are independent of S, we can take the limit as 5 — > 
in last equation to obtain 



E, 



J l K (X u )U Xu (A £ )dL* u = E X , 8 [ l K (X u ) Cl U Xu (A")dL v 



A standard argument involving the monotone class theorem, shows that this last 
equation is not only valid for Ik , but also for all bounded, measurable functions 
/ : dD —> K. Since for each e > 0, the function x h- > H X (A £ ) is bounded away 
from zero (see [4]), we can take f(x) = e^ at H x (A £ )^ 1 to get 



e~ at dL* u =E x , s e- at Cl dL u , (4.13) 
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for any positive a, and arbitrary v < t. We can extend (4.13) in the following 
way: for N € N and T > 0, define simple functions by 



N 



k=0 



It follows from (4.13) that E XjS J °° fjy t T(u)dL* — E X)S J* °° fN,T{v)c\dL u . It is 
clear that /jv,t(w) increases to e _QU l[Q !T ] (tt) for all T > 0, and so, by the 
monotone convergence theorem, we obtain 



E x . / e~ au dL* u = E XtS / e~ au Cl dL u 



that is, L* and L have the same a-potential functions. Since both L* and L are 
continuous, it follows by [3], Chapter 4, Theorem 2.13, that c\L = L* a.s. This 
shows that (dL t , c±H x ) is an exit system, an the theorem is proved. 

It remains to show that equation (4.12) holds. It might seem that the exit 
formula (3.14), and equation (4.12) just differ on the upper limit of the integral. 
In the former, the upper limit is a fixed time T whereas in the latter, the upper 
limit is the random time t j . The advantage of having r J as an upper limit is 
that at most one excursion that reaches level e away from dD starts before time 
t j . We will drop the super index j from t ] for notational simplicity. 

Call I z to the left hand side of (4.12). The process l{ r<u }, and thus 1{ U < T }, 
are predictable. Since at time r the process X cannot be on the boundary, we 
see by the exit formula (3.14) applied with V u = l{ u < T }lif J . 

/>oo 

h=^z Mu<r} 1k 3 (X u )U Xu {A e )dL* u 
Jo 

/•oo 

J° u<oo 

= E 2 hu<T}lK j ( X u)~LA*(eu)=® z J2lK j (X u )t A e(e u ), 

which roughly says that the exit formula is also valid if we change T > for the 
stopping time r. 

Since only one excursion in A e can happen before time r, we see that I z is 
equal to the probability that, after escaping Cj, an excursion reaches distance 
e away from dD. Intuitively, a reflected Brownian motion Y starting at z accu- 
mulate roughly the same local time as X before exiting Cj. Since the cylinder 
Cj is very thin, both Y and X are likely to exit the cylinder through Cj and 
be away from each other no more than the amount of local time accumulated 
up to time r. Since both X and Y are the same Brownian motion inside of 
D, the probability that, after exiting Cj, an excursion of X reaches distance e 
from dD is roughly the same as if such probability is computed with respect 
to Y . This idea yields that I z can be estimated by using Y instead of X . For 
reflected Brownian motion, it is known that L* can be chosen to be the local 
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time from its Skorohod decomposition, that is, Y satisfies (4.12). So, X should 
satisfy (4.12) as well. We will formalize this idea next. 
Dchnc A = {3 u < t, X u € Kj, e u £ A £ }. Then, 

^J2lKA x unA4eu)=V z (A)=F z (A,X T e Cj) (i + pj^ecf) ) ' 

Call D £ to the set of points in D at distance at least e from dD. Let Toe be 
the hitting time of D £ and the hitting time of Cj. By the strong Markov 
property, 

P z (k,X r £ Cj) _ P. (Px T (Ti>. < T dD )l dcACJ (X T )) 
P Z (A, X r i Cj) ~ Pz U {Tdc ~ Tdo)lcT {Xt[ 



The process B t = B (T^e — t) is also a Brownian motion by the strong Markov 
property and independence of increments. Therefore, 

p z (A, x T i cj) sn VxeD . v x (b TC] e dc, \ cj) u y{dCj \ cj) 

^ : S: C E 



¥ z (A,X T eCj)~ M xeD , P x (B Tc . G Cj) ~ 6 ujV(Cj) ' 

where uj v represents the harmonic measure of D \ Cj, for an arbitrary point 
y € D £ , and C e > only depends on e. The surface area of the side of the 
cylinder Cj is about VS times less the surface area of the top of Cj, so we 
conclude that 

h = (1 + 0(S))¥, (Px T (lz> < T 9d )1 c t(X t )^ . (4.14) 

Next we introduce a reflected Brownian motion Y starting from z, driven by 
the same Brownian motion B that drives X. We proceed to do some estimate 
to compare X and Y. 

We claim that P Z (L T > 2r$\/~S) = O(S). Otherwise, there is a sequence 5 n I 
and p > such that P Z (L T > 2rs n y/6^) > p. But any function u> e C[0, 1) with 
L t (uj) > 2rsVS satisfies that B^(lo) < for t < Tg(uj). This only occurs with 
probability 0(5), which contradicts the existence of the sequence 5 n - Since 

\\X t -Y t \\ 2 - \X? y/f < 5 2 (L T + L Y T f + HIlL L 2 T , 

the process Y leaves the cylinder Cj through the side only when X t gets at 
distance (S 2 + r$^/5)0(l) from the side of Cj. This event has probability of 
order ^/S by a harmonic measure argument similar to that that lead to (4.14). 
Therefore, the event Y u g Cj for u < r has probability 1 — 0(8). Actually, the 
event holds true if we request that \z — Xj\ < rs — 2rsVS, as long as L T < 2rsVS. 
Thus, since for y E Cj Pi dD we have 1 > n(y)e n > 1 — S 2 , 

X?-Y t n =( n(X u )e n dL u - f n(Y u )e n dL Y = (1 + 0(5))[L t — Lj], 
Jo Jo 
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for all t < t, on the set |i T < 2rg\f5^- 

It is clear that dist(X T , 3D) is comparable to rsVS, and that dist(Y T ,dD) 
is comparable to Y?, since Y™ does not leave the cylinder Cj on the event 

Irs We will argue that \X™ — Y™\ < Ar$5. Indeed, assume that 
X™ > 2r§8. Since no piece of the boundary is above level 2rs8, no local time L 
is accumulated and so X n — Y n does no increase. Similarly Y n — X n does not 
increase when Y n > 2r$5. But if both X n and Y n are less than 2r$S, then it is 
clear than \X n — Y n \ < 4r$5, since the piece of boundary Cj n <9D is between 
levels — 2rsS and 2r$S. Thus, 



^- - 1 




yn _ X n 






X? 



and so, for some positive constants as, b$, depending on 6, 

dist(Y T ,&D) 



as < 



dist(X T ,dD) 



< b&. 



Thus, by the boundary Harnack principle, there are positive constants cs, Cs 
such that 

?X T {Tjje < T dD ) 



cs < 



<C S . 



Py r (T D . < T dD ) 

Moreover, since \\X T — Y T \\ — >• as S — ¥ 0, Lemma 1 in [5] shows that both cs 
and Cs can be chosen arbitrarily close to 1 as 6 — > 0. Using this fact in (4.14) 
we obtain 

I z = (1 + 0(5))P Z (Py T (T De<TaD )l C3 (y r )) (4.15) 
+ (1 + 0{S))¥, (f Xt {T d . < T 9d )1{ Lt 

>2r s Vs} ) ■ 

Since dist(X T ,dD) is about 2rsV~5, standard estimates for Brownian motion 
show that Fx t {Td s < Tqo) is comparable to e~ 1 rs\ /f 5, so 

= 0(5)E Z (L T ). 

By using this last equation in (4.15), we are allowed to write the following 
equation, by the same arguments that led to (4.14), 

I z = (1 + 0(S))F Z (¥ Y ^T D , <TgD )l Cj (Y T )) + 0(S)E x (L T ) 
= (l + 0(5))E z (^l c .(F u )l Ae (e^)) +0(S)E x (L T ), 



where e Y denotes excursions of the reflected Brownian motion Y. Theorem 7.2 
in [4] says that (dLj , ciH y ) is an exit system for Y for some constant ci > 
independent of y, therefore 



(l + 0(S))E z 



f l Cj (Y u ) Cl nY u (A e )dLl 
Jo 



+ 0{S)E z (L T ). (4.16) 
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We will use the Harnack boundary principle for the harmonic function x 
Pf (A £ ), which vanishes on Cj n 3D. Since 



By the Harnack boundary principle, and Lemma 1 in [5], we obtain that the 
right hand side above is of order 1 + 0(5). It follows that 



E 



[ T l Cj (Y u )U Yu (A s )dLl = (l + 0(6))U Xj (A e )E z (L^). (4.17) 
Jo 



The same argument allows us to write 
E 



/ l Cj (X u )n Xu (A e )dL u = (l + 0(6))H Xj (A e )E z (L T ). 
Jo 



(4.18) 



By the previous estimates and the optional sampling theorem applied to 
E,(L T ) _ {i + om E z (X?-B? + 6*L T ) 



(1 + 0(6)) 1 + 



E,(r r "-B«) 
e,(x? - + 6 2 L T ) 



But, 



e z (x™-f™ + £ 2 £ t ) 



: (y» -x? + x») 



< 



E z (y™ - x? + x? 

Ar s 6 + 6 2 E Z {L T ) 



-4rsS + 2rs\/6 



0(6), 



which shows that E Z (L T ) = (1 + 0(5))E Z (L^). Using this fact, equations (4.16), 
(4.17), and (4.18), we obtain (4.12), as desired, and the theorem is proved. □ 



4-3. The stationary distribution 

Lemma 3.21. Let D, f and g be as above. Let r,T > 0, and z £ D. Assume 
that Al holds. Then, for any (x ,s ) £ D x W there is uj £ C([0, T]; R") with 
bounded variation such that there is a unique (x, s) £ C([0,T]; DxM?) satisfying 
(x(0),s(0)) = (x ,s ), (x(T),s(T)) = (z,0), and for t £ [0,T]. 

x(t) = xq + uj(t) + / 7(x(u), s(u))dl(u), 
Jo 

k(t) = fco + / g(x(u)) — s(u) dl(u), 
Jo 

Here, l(-) is a continuous and increasing function, that only increases when 
x(t) £ 3D, that is l(t) = J * l dD {x(u))dl(u). 
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Proof. In virtue of Theorem 4.1 in [16], we would expect uniqueness to hold in 
the bounded variation case. The issue to apply such theorem directly is that 
the reflection vector 7 depends on the value of fc(-), but it is easy to adapt the 
proof of uniqueness in such theorem to our setting, just like we did in the last 
step of Lemma 3.15. 

Next we construct a function ZJ £ C([0,T];R") with bounded variation, and 
a solution (x, s) of the system above. Consider the uniform partition < a\ < 
h < a 2 < ■■■ < b p+1 < T of [0,T]. 

To construct ZJ and the associated solution, set uJ(0) = and for t £ (0, ai], 
let ZJ(t) be defined as any fixed continuous function with bounded variation, 
such that Xo + u(t) £ D, and xq + oj{a{) = x\. It is clear that any solution (x, s) 
has to satisfy x(t) = xq + ZJ(t), s(t) = s Q , and l(t) = up to time a\. Next we 
want to keep x{t) at x\ from a\ to b\. In view of (3.17), for t £ [01,61] we set 
= exp [l(t) - l(ax)} = e 1 ^ and so 



By setting l(t) = + r/i(t — a{) for t £ [a\, 61], where 771 is to be determined, 
we obtain that both I and s are continuous. All this implies that we need to 



Uniqueness in [ai,£>i] follows directly form the fact that the equation above 
defines a continuous function with bounded variation. Thus, the functions (x, s) 
defined above correspond to the unique solutions to the Skorokhod problem for 
ZJ in [a%, bi]. 

Next, we want to keep s(t) constant in [61,02], while we move x(t) from 
xi to X2- To this end, pick a curve (1 : [61,02] — > D such that Ci(^i) = x i- 
Ci( a 2) = %2 and (i(t) £ D for other values of t. Set l(t) = l(bi) for t £ [61,02], 
and ZJ(t) = ((t) — x\ + ZU(6i). It is clear that the only solution with bounded 
variation in this interval is x(t) = £i(i) and s(t) = s(6i). 

We iterate this process by keeping x(t) at Xj in [aj,bj), and by defining 
l(t) = I (bj-i ) + r/j (t — a,j ) , yj (t) = exp[77j (t — aj )] in that interval. This way, the 
function s(t) must satisfy 



s(t) = yiityho + gixMty 1 I e l ^dl(u) 



= y^tyKso+gixMi)- 1 [yi{t) - 1] . 



define 





(4.19) 



j 



+ X] 9( x m)[y m (b m ) ■ ■ ■ yj-i(bj-i)yj(t)] 1 [y m (b m ) - 1] 



for t £ [a,j, bj]. The calculation leading to such equation, though tedious, is 
straight-forward to carry out by splitting the integral in [0,t] into integral in 
the sets [aj 7 bj] and [6j,6 J+ i], and using our definition of l(u) in each of those 
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intervals. In the interval define u(t) by 

uj(t) = u{a 3 ) - rjj I j(xj,s(u))ds. 

Once again the unique solution in this interval for this uj is (xj, s(t)). 

From bj to ay+i we find a curve Q going from Xj to Xj+i through D, and set 
tij(t) = (j(t) — Xj+uJ(bj) and l(t) = (bj). The unique solution is then (£/(*)> s(bj)). 
This procedure can be also done so that x(T) = z. 

It remains to show that we can choose the values of rjx, . . . , rj p +i > such 
that s(T) = 0. At time T we find that 

P+i p+i p+i 

s(T) = s Y[ DmibmT 1 + g( x ™)[y™{b m ) - 1] II ViihY 1 , 

m—1 m—1 i—m 

and to obtain s(T) =0we need 

p+1 m—1 

-so = ^2 g( x m)[y m (b m ) - 1] JJ Viih). 

m—1 i — 1 

By Al there are non-negative Ai, . . . , A^+i such that — sq = Ylm=i d( x m)^m, 
so we need to choose the numbers r\ m so that y m (b m ) satisfies A m = [y m {bm) — 
1] nr^i 1 Ui(bi)- This is easily achieved by an inductive procedure, and the lemma 
is proved. □ 
Theorem 3.24. Let (X,S) be spinning Brownian motion solving (1.1) with 
So = 0. Then, for every T > 0, there is an open set U C H^, and c > such 
that for every open B C U it holds that F x , s {St € B) > c m p (B). 

Proof. In Theorem 3.17 we have obtained an exit system for Z t — (X t ,St) 
representing excursions from dD xR p . As St does not change within an excursion 
of Z t away from dD xW, we can think of the excursion law H I S as a measure 
representing paths of X only, and thus we will drop the subindex s. 

From now on we closely follow part of the proof of Theorem 6.1 in [2]. To sim- 
plify the notation, call Z — (X, S), and we use the standard nomenclature 
for first hitting time of a set U by the process X, and a t = inf {s > : L s > t} 
for the right inverse of local time. 

We proceed to describe an exit system for a different, though related, process 
X' . Let zo G D and r > be arbitrary but fixed, so that B r (xo) Q D and 
set U = B r (zo). Let X' = X' u be the process X conditioned by the event 
\Tjj > cti}. It follows from Proposition 3.22 and the strong Markov property 
that for any starting point in D, the probability of {Tjj > ci} is greater than 
zero. It is easy to see that (A t ', L t ) is a time homogeneous Markov process under 
P ZOjSo in (y^t^o- To be consistent with the notation, we will write (X' t ,L' t ) 
instead of (X[,L t ). 

The arguments in step 2 in the proof of Theorem 6.1 in [2] can be followed 
without any essential modification to use the exit system (dL t , H x ) to construct 
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a new exit system, for the process (X' t u ,L' t ), and use it to show how X t up to 
time can be recovered from X' t u by adding a las excursion that hits U. In 
particular, there is a measurable map T : C x 23(£>) — > G such that T(X' U , U) 
and X have the same distribution up to time . 

Let Uj = B r (zj) for j = 1, . . . ,p + 1, where zj G D are chosen so that the 
sets [/j are pairwise disjoint and their union is a subset of D. 

Let X 1 be the process X' U2 satisfying that X$ is uniform in U±. Define Y 1 = 
T(X 1 , XJq). The process Y 1 is an sBm starting with uniform distribution in U\, 
observed until the first hit of U2, at time T\ = inf {t > : Y^ G t^}- Similarly, 
for j = 2, . . .,p we define X 3 to be the process X' Uj+1 starting with uniform 

distribution in Uj, and set Y j = T(X j , U j+1 ) and Tj = inf |t > : if G £^+i}- 

It should be clear that the processes can be chosen to be pairwise independent 
and they have the distribution of a sBm starting with uniform distribution in 
Uj observed up to the first hitting time of Uj+\. 

In the exact same way as in Step 3 of the proof of Theorem 6.1 [2], we can use 
the processes {Y^} to construct a spinning Brownian motion X* in D, starting 
with uniform distribution in U\, and such that, conditional on {Xj,t > 0}, 
j = 1, ...,p + 1, there is eg > such that with probability at least eg we 
have that X* is a time-shifted path of X\ on some appropriate interval, for all 
3 = 1,... ,p. 

To conclude this proof, we show that with a positive probability, the process 
S can have "almost" independent and "almost" linearly independent increments 
over disjoint intervals of time. This is used to show that a conditioned version 
of S has a density, or equivalently, that S has a component with a density, in 
an appropriate open set. 

Let X\, . . . , Xp+\ be points in dD satisfying assumption Al. Since the matrix 
• • • |g(:Ep+i)] has rank p, it is possible to eliminate a column from it and 
still have a matrix with rank p. It follows, withouth loss of generality, that 
the vectors g(x±), . . . ,g(x p ) are linearly independent. For 1 < j < p, let Cj = 
{z £ K" : Z(g(xj), z) < 5q}, for some 5q > so small that for any Zj € Cj, 
j = 1, . . . ,p, the vectors {zj} are still linearly independent. Let 5i > be so 
small that for every j = l,...,p, and any x € dD n B$ 1 (x 3 ■), we have g(x) G Cj. 

Let V the local time of X^ on dD and a{ — inf {s > : U s > t\. It is not 
hard to see that for some P2 > 0, the probability that for every j = l,...,p we 
have X\ ^ dD \ Bg 1 (xj), for t £ [0, al], is greatr than p 2 . Let 

Rj = sup |t < Tj : : Y t j G 8d\ and Qj = L j R . . 

Consider the event -F* containing all trajectories such that for j = 1, . . . ,p, 
we have X\ £ dD \ Bs x (yj) for t G [0, a{] and Rj < a\. The construction of the 
process X' in [2] then shows easily that for some constant C3 > the inequality 
P 20 (^) >P2(l-e- c *)P holds. 
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Let Ej = exp ( L\ ) , and define the following collection of random variables 



Notice that if holds, then S 3 (t 3l . . . ,t p ) £ Cj for all tj,...,t p € (0, 1] and 
j = 1, . . . ,p. For any < a k < b k < Qfc for fc = 1, . . . ,p, define the random set 



A([ai,6i], . . . , [flip, ftp]) = < ^2 S j (tj, . . . ,t d ) : t k G [a fe ,6 fc ] 

It is not difficult to estimate the d-dimensional volume of A([ai, &i], . . . , [a p , b p \) 
by using the definition of Cj's. First, by continuity, it follows that under F+ 
there is a positive constants q such that (1 — qSo)e tj < E d t < (l + q5o)e tj for all 

1 < j < P- By definition of the set Cj , it follows that for some positive constant 



(1 - pS )g( Xj ) - 1] < / J g(X^)Erdu < (1 + /Wofe) [e*' - l] . 

« o 

Define a function t?(ti, . . . , tj) by 

j/=^J e *=i [e*> - 1] . 

i=i 

It follows from the inequalities in this paragraph that for some constant rj > 0, 
independent of 5o, we have 

p 

(1 - TjSoMh, . . . ,t d ) <J2 s Htj, ■ ■ ■ ,td) < (l+r)5 )v(ti,...,t d ). 

3=1 

If is small enough, we have that there exist a constant C3 independent of 
a k , 6^ such that 



-1 < m d (A([ai,bi],...,[a p ,6p])) < 
' :; " m d {v(t u ...,t p ) : t fe e [a k ,b k ]} ~ ' ' 



To compute the volume of the set in the denominator, we will calculate the 
Jacobian of v(-) in Lemma 4.1. We obtain that det Dv = C exp (— X)fc=i kt k ), 
which, as t k € [0, 1], readily yields that the d-dimensional volume of the random 
set A([oi, 61], . . . , [a p , bp]) is bounded above by 04(61 — ax) •• ■ (b p — a p ) and below 
by 05(61 - 01) • • • (b p - a p ). 

Let us consider the processes X* defined above, conditioned on the sigma 
field 

S = a({S j (t j ,...,t d ), t k e[0,l}} d j=h 
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The construction of X' in [2] makes the random variable Qj the time component 
of a time-excursion Poisson random variable with variable intensity given by 



where a' is the right inverse of the local time of X' and A is a fixed set such that 
C5 < H X (A) < cq for positive constants cs,C6 that are independent of x G dD. 
It follows that Qj has a density with respect to the Lebesgue measure in [0, 1] 
that is bounded below. In view of our remarks on the volume of A, it follows 
that conditional on S, the vector 



has a density with respect to the p-dimensional Lebesgue measure that is bounded 
below by C7 > on the open set U = A((0, 1), . . . , (0, 1)). We can now remove 
the conditioning on and conclude that S(Qi, . . . , Qd) has a component with a 
density with respect to p-dimensional Lebesgue measure that is bounded below 



on U. 

Define E* = exp(L*) and S* t = E* t ~ 1 J* g(X*)E* dL* s and T* = £^7}, 



where L* is the boundary local time for spinning Brownian motion X*. Using 
conditioning on F+, we see that the distribution of has a component with 
density greater than eg on U. 

The previous argument can be modified to show that for any fixed to > 0, 
the random variable S* o , 2 has a component with a strictly positive density with 
respect to d-dimensional Lebesgue measure on a non-empty set U, which proves 
the theorem. All that we need to do is, for small e > 0, find times tj > such 
that Tj e (tj — e, tj + e), with uniformly (in j) positive probability p e , and then 
further condition the stitched process X* to satisfy Tj e (tj — e,tj + e). Set 
** = This way, T* = J2j=i Tj € (i* — ed, t* + ed) and since G Ud, 

and Ud is away from the boundary, we can condition on A* to not to hit dD in 
[U - ed, t* + Ed] and thus have = . We then choose t = . □ 

Lemma 4.1. Let v(-) be the function defined in the last step of the proof of 
Theorem 3.24- There exists a constant C / 0, depending only on the vectors 
g(xk) such that 




S(Qi, ...,Q d ) = S 1 (Q ll ...,Q d ) + S 2 (Q 2 , ...,Q d ) + --- + S d (Q d ) 



d 




for all t!,...,t d e (0,1]. 

Proof. A stright forward calculation shows that for i — 1, . . . , d: 
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Let Xj i the coefficient of g(xj) in the expansion of ^f- above. Because the 
vectors g(xk) are linearly independent, these numbers are well defined. We spe- 
cially remark that Xj^ = for j > i. Let Tg the dxd matrix whose j— th column 
is g(xj) and let A be the matrix whose component is Xj^. The calculation 
above then simply says that Dv = TgA. Therefore, as A is triangular 



The double sum on the right hand side equals X^fe=i ^ifc, by Fubini's theorem, 
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